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ABSTRACT 
This  thesis  is  the  result  of  literary  research  on  the  subject  of 
sensitivity  in  linear  feedback  control  systems.   It  is  a  synopsis  of 
information  obtained  from  various  technical  publications  and  is  designed 
to  give  the  reader  information  on  the  theory  of  sensitivity  and  its 
application  in  design  problems.   Included  are 

(1)  Definition  of  sensitivity  function  and  root  sensitivity 
for  incremental  variations  of  single  parameter  and  their 
design  application. 

(2)  Definition  of  sensitivity  function  applied  to  the  large 
parameter  variations  and  its  design  application. 

(3)  Special  topics  such  as  specification  for  sensitivity 
function,  sensitivity  integrals,  and  introduction  to 
multi-parameter  sensitivity. 
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INTRODUCTION 

One  of  the  foremost  properties  of  feedback  is  its  ability  to 
reduce  the  sensitivity  of  a  system  to  variations  of  the  system  para- 
meters.  Suppose  there  is  an  element  in  the  system  which  must  be  used 
because  of  its  special  properties.   However  the  element  is  sensitive 
to  environmental  conditions  (such  as  temperate  re  and  pressure),  or  is 
subjected  to  aging,  or  has  wide  manufacturing  tolerances  such  that  when 
it  is  replaced  the  parameters  of  the  new  element  may  be  markedly  diff- 
erent from  those  of  the  old  element.  Whatever  the  reason  may  be, 
suppose  the  values  of  the  element  parameters  vary  significantly.   It 
is  desired  that  despite  these  variations  certain  system  properties, 
such  as  the  input-output  transfer  function  or  output  impedance  should 
remain  substantially  constant.  The  troublesome  element  could  perhaps 
be  replaced  by  a  better  one,  but  this  may  be  very  expensive  or  imposs- 
ible.  In  any  case  we  would  like  at  least  to  ave  an  alternative  to 
replacing  the  element.   In  most  cases  the  troublesome  elements  are  the 
active  elements,  vac  m  tubes,  transistors,  energy  conversion  devices 
such  as  motors  and  generators.   :.■  t  tbey  may  also  be  passive  elements 
or  transducers.   In  general  the  system  assumed  in  an  analytical  study 
will  never  precisely  match  the  actual  physical  system;  it  is  therefore 
important  to  know  the  effects  of  possible  variations. 

A  general  knowledge  of  the  effects  of  parameter  changes  can  be 
used  as  a  guide  to  system  modifications  which  would  improve  the  over- 
all performance.   For  all  these  reasons  measures  of  closed  loop  system 
sensitivity  to  open  loop  system  parameter  variations  are  an  integral 
part  of  the  analysis  problem. 


Notions  about  system  sensitivity  were  in  the  forefront  when  the 
feedback  concept  was  initially  developed.  This  was  natural,  even 
unavoidable,  since  feedback  systems  possess  the  fundamental  physical 
property  that  the  effects  of  variation  in  the  forward  loop,  whether  they 
are  taken  as  changes  in  open  loop  transfer  function  G(s),  or  as  depart- 
ures from  strict  linearity  or  from  freedom  from  extraneous  noise,  are 
reduced  by  the  factor  1  +  G  in  comparison  with  the  effects  which  would 
be  observed  in  a  non  feedback  system.  Accordingly,  sensitivity  measures 
were  Indispensable  to  any  rational  discussion  of  feedback  systems  and 
a  useful,  classical  definition  of  sensitivity  was  made.  Except  for  a 
minor  modification  the  definition  of  sensitivity,  given  by  Bode,  remained 
unchanged  for  over  a  decade.  Perhaps  this  static  nature  surrounded  by 
dynamic  growth  in  most  other  areas  of  feedback  systems  engineering, 
made  the  concept  fade  with  little  exposure  to  sensitivity  concepts. 
In  recent  years  classic  sensitivity  has  become  a  more  popular  subject 
in  automatic  control  and  also  in  network  synthesis.  Finally,  the 
emphasis  on  pole-zero  specifications  for  system  characteristics  gave 
rise  to  new  concepts  of  sensitivity,  with  associated  new  measures. 
These  measures  called  here  "gain,"  "open-loop  pole,"  and  "open  loop 
zero"  sensitivities,  were  evolved  to  account  for  changes  in  the  posi- 
tion of  closed  loop  poles  due  to  shifts  or  changes  in  open  loop  gain, 
poles  or  zero. 

This  thesis  is  the  result  of  a  literary  research  on  the  subject 
of  sensitivity  in  linear  feedback  control  systems.   In  section  II,  the 
system  sensitivity  and  root  sensitivity,  for  a  system  subjected  to 
incremental  parameter  variation,  are  discussed.  Based  on  the  meaning 
of  open  loop  pole  and  zero  sensitivity,  a  design  technique  for 
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compensating  the  system  is  also  introduced.   In  Section  III,  the  large 
parameter  variation  and  a  new  definition  of  sensitivity  is  presented. 
It  is  shown  that  when  the  leakage  transmission  between  input  end  output 
is  negligible,  the  new  definition  and  the  classical  definition  of  the 
system  sensitivity  are  identical.   In  Section  III,  an  attempt  is  made 
to  utilise  Nichols  Charts  for  evaluation  of  the  sensitivity  at  a  part- 
icular frequency  range.   A  design  philosophy  and  technique  is  also  dis- 
cussed.  Section  IV  discusses  some  special  topics  which  are  of  parti- 
cular interest  today  and  could  be  selected  as  interesting  topics  for 
further  investigation. 
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II 

SENSITIVITY  AND  INCREMENTAL  PARAMETER  CHANGES 
II-l  -  System  Sensitivity  and  Root  Sensitivity 

Several  definitions  of  sensitivity  have  been  offered  in  the  past. 
Bode  defined  the  sensitivity  of  the  overall  transfer  function  T  to  the 
plant  parameter  k  as 


JlL 


s 


k     \jU1  \1 


Horowitz  used  the  inverse  of  Bode's  definition, 


gT  _       hJUT 


*T 


k  &/*k  frk 

k 

T 
Defined  one  way  or  another,  S.  is  generally  known  as  "classical  sensi- 
tivity" or  "system  sensitivity,"  since  it  involves  the  system  transfer 
function  T.  Several  other  definitions  have  also  been  suggested,  and 
each  one  of  these  definitions  concentrates  on  some  particular  engineer- 
ing or  mathematical  simplicity.  A.  Schulke  suggests  since  sensitivity 
relates  changes  in  the  transfer  function  of  a  system  with  respect  to 
changes  in  parameter  of  the  system,  then  one  can  define, 

i m  _i_  m       fr  T 

sensitivity       S   "   ^x 
and  in  general  T  is  a  complex  and  may  be  expressed  as 

T(x)  -  A(x)  +  jB(x) 
where  A(x)  is  the  system  attenuation  function  and  B(x)  is  the  system 
phase  function,  and  from  that, 


-L.  -       U       .  _Jl_ 

S      IjUx  >>x 


x 


j  A(x)  +  j  B(x)j 


*A(x)  oB(x) 

*x    +  JX    fcx 
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Then  assuming  that  x  is  real ,  it  is  apparent  that  the  real  part  of 

;r  is  directly  related  to  the  system  attenuation  and  that  the  imaginary 
s 

part  is  directly  related  to  the  system  phase  function.  The  use  of 
S*  which  Schuike  calls  it  "sensitiveness"  permits  comparison  of  both 
attenuation  and  phase  functions  relative  to  changes  in  a  common  para- 
meter 

The  most  adopted  definition  of  system  sensitivity  these  days  is 
the  one  suggested  by  Horowitz,  i.  e. 

and  we  adopt  it  through  the  coming  discussions. 

In  recent  years  with  the  increased  utilisation  of  the  pole-zero 
approach,  it  has  become  increasingly  important  to  examine  variations 
in  position  of  the  poles  and  zeros  of  the  network  function  due  to 
changes  in  network  parameters.  This  is  known  as  the  "root  sensitivity." 
Formal  definitions  of  root  sensitivity  vary  from  author  to  author. 
.Some  define  it  as, 

s .       *s .  ^8 


;  i   »        3     ,   x       J 
X  OX  0  X 


where  (-s.)  is  the  root  of  the  characteristic  equation  and  x  a  para- 
meter which  could  be  the  gain  K,  or  a  pole  or  a  zero.  Huang  uses  the 
definition 


sj 

s  J 

- 

b8J 

- 

X 

ox 

X 

&* 

M  Fuer  ad  Stapleford  prefer  different  definitions  for  sensitivity  to 

SJ  Sj  FJ 

gain,  S_  ,  and  sensitivity  t     poles,  S  ',  or  sensitivity  to  zeros,  S_  . 


Again  each  one  of  these  definitions  leads  to  a  particular  engineering 
or  mathematical  simplicity  utilized  by  that  particular  author. 
II-2  -  Sensitivity  Relation  in  Control  Feedback  Systems 

It  has  been  noted  by  Bode  that  for  any  linear  system  the  transfer 
function  T(s,K),  relating  a  response  to  an  input  function  can  be  in 
the  bilinear  form  of, 

Tfs  K)   -    A(s)  +  KB<s> 
ll*>*>  C(s)  +  KD(s) 

and  T  m       )faT         .   -^L 


K 


Both  numerator  and  denominator  of  T  are  of  the  form, 

H(s,K)  -  q(s)  +  Kp(s) 
where  K  is  a  variable  parameter  and  p(s)  and  q(s)  are  polynomials.  The 
roots  of  H(s,K)  »  0,  which  characterise  H,  depend  on  the  value  of  K. 
H(s,K)  «  0  constitutes  an  implicit  relationship  between  s  and  K,  and 
as  K  takes  on  all  positive  values,  a  number  of  curves  in  the  S -plane 
are  obtained.  These  curves  are  the  loci  of  the  roots  of  H  called  the 
"root  loci." 

The  root  locus  was  introduced  by  Evans  in  1948.  Re  used  it  main- 
ly for  the  investigation  of  control  systems  with  a  transfer  function 
T,  in  the  form  of, 


T   ■ 


KP(S) 


1  +  G         q(s)  -I-  Kp(s) 
where  G  -  Kp/q  is  called  the  forward  gain  in  a  unity  feedback  con- 
trol system.   For  such  systems,  only  the  denominator  is  of  interest. 
As  the  result,  there  arose  the  generally  accepted  terminology  in  which 
the  zeros  of  p  and  q  were  called  open  loop  "poles"  and  "seros," 
respectively  and  those  of  1  +  G  ■  0  are  called  "closed  loop  poles." 
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II-3 


-  Relationship  Between  System  Sensitivity  and  Pole  Zero 
Sensitivities 

The  transfer  function  T(s)  can  generally  be  written  in  the  form 


of 


n 


T(s) 


6irr(s  +z.) 

X 

m 

TT<8  +  8.) 

1      J 


(2  1) 


and  let  Z  ,  s  and  g  be  functions  of  some  parameter  x  (such  as  temper- 


'.!■ 


ature,   gain  of  some  element,   etc.).     Then 

n  m 

JUT    -    /*g    +   Xuts  +   Z  )  -  2M»  +  s4) 

1  J  1 

and  differentiating  above  we  get 
>  dx 


dT 


g  ax    1  >x   s  +  Z.    1  dx   S  +  8. 


g 


dividing  both  sides  by  dx 

x 


dT 
T 

dx 
x 


i£ 


*Z 


*2      i    bx  s  +  Z 


m 


2   is. 
x   J 


(2-2) 


1   ix  s  +  s. 


Equation  (2-2)  gives  the  relationship  between  the  classical  definition 
of  se  sitivlty  and  the  closed  loop  pole  and  zero  sensitivity,  where 
the  closed  loop  pole  and  zero  sensitivities  are  defined  as 

s  v 


■  x 


*x 


e«*  •  Zx 

X        *x 


(2-3) 


(2-4) 


This  definition  is  not  necessarily  the  most  useful  definition  at  all 
situations. 

Pole  Sensitivity  (closed  loop  pole): 
T> e  root  locus  is  defined  by: 
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H(tf,R)      -      q(s)     +     Kp(s)      -      0 


(2-5) 


hence- 


TT<8+  p  ) 
1»1  J 

k   -        <*(«)      " 


P(8) 


m 


TT(8  +  Z  ) 
1-1 


(2-6) 


and 


rs 


m 


jum    -    ZMs  +  p«>  -  £A<8  + z.)  + 

j=l  J       i=l 


T^J 


(2-7) 


and 


dR 

K 


44 


n  -  m  l 


Also  since 


(2-8) 


then 


and 


K     -     - 

P(8) 

dR 
ds 

q'p   -  P*q 
P2 

dR 

ds       " 

q'p  -  p 

2 
P 

q 

R 


a_  -  £_ 
q      p 


(2-9) 


(2-10) 


For  simplicity  of  notation  Equations   (2-6)  and   (2-8)   can  be  written  as, 

.A, 


K     -     .TT(i  +  Si)*i 


*K  A 

1     "     dS 


i     s  +  s 


(2-11) 


(2-12) 


which  has  to  be  evaluated  at  s  ■  -s..   The  roots  of  H  ■  0  are 

J 

denoted  as  (-s.)  to  distinguish  them  from  the  variable  s  which  may  of 
course  assume  any  value.   In  other  words  -s ,  is  a   point  on  the  root  locus, 
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And 


if  -  s .is  a  root  of  q(s) 


if  ~  s  is  a  root  of  p(s) 


-4  -  -l 

From  Equations  (2-10)  and  (2-12)  we  can  write, 

8 


J 


ds 


K      i 


8  +  8 


j£3- 


qp'  -  pq 


S»8 


i»  ju 


Kp' 


(2-13) 


where  S   is  the  closed  loop  pole  sensitivity  with  respect  to  K  and 

sometimes  referred  to  as  just  pole  sensitivity  with  respect  to  K. 

Furthermore  the  residue  of  a  function  A(s)  at  pole  s»a  is  defined  as 

B(s) 


Residue 
at  s=>  a 


B(s) 


lim   (s  -  a)  A(s) 
s-»a      B(s) 

lim   (s  -a)  A'(s)  +  A(s) 
s-a  "      B'(s) 


-  lim 


MsL 


s-»a   B'(s) 

which  is  true  for  simple  poles.  Noting  that  feg   »  q'  +  Kp'  and 

Si" 
applying  above  residue  theorem  we  get 


J 


Pj 


K 


kl 


q'  +  KP' 


where  Q  is  the  res id  e  of 


T(s) 


Kp(s) 


q(s)  +  Kp(s) 


J&, 


>8j 


(2-14) 


at  point  8  ■  -s  .  The  same  result  can  be  obtained  directly  from  (2-5) 


ds 

dK 
K 


-   -  K 


Ml 


(2-15) 
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and  since  H  ■  q  +  Kp,  then 


'K 


-  K   fcK     -   -      Kp 

-  -  Kp 

3H           dq  +  K  dp 

fts          ds     ds 

s.-8j 

SRJ  -   -Kp/ifq*  +  Kp*)  is  quite 

convenient  sine 

-o, 


S--6 


The  relatii 

requires  knowledge  of  the  "open  loop"  transfer  function  and  need  not 

be  in  factored  form.  The  evaluation  of  the  pole  sensitivity  S_J  by 

the  residue  of  the  closed  loop  transfer  function  will  call  for  the 

complete  determination  of  the  closed  loop  transfer  function. 

II-4  -  Slope  of  Root  Locus 

Slope  of  the  root  locus  is  the  tangent  of  the  argument  of  & 


and  since, 

c8J   . 


4s 


1 


AK 
K 


s 


Then  As 


J 


J 


K    K 


j 

(2-16) 
(2-17) 


Since  AK  is  a  positive  number,  then 

Ab.     -  /sj 


where /As.   represents  arg  As. 


dx 


tan  /Ls  .     *»     ten  /  s„ 


So 

tsi 


w 


KeFvJ 


8  . 

K" 


(2-18) 


(2-19) 


Vk 


Kj 


1  (2- 

•8)  we  have, 

1 

I 

m 

2   * 

m 

-  2  1 

S»-8 

i               - 

Be 

2       * 

[A  4  +  pl 

m     x 

8«-8 

(2-20) 
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Input 


-**©- 


Output 


Fig.  2-1:   Unity  Feedback  System. 


8  *  x  +  jy 


Fig.  2-2:   Root  Locus  Slope. 
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If  we  Write 

.  +  »j  ■  i,  •** 

and  since 


I     /        1  \        I  1  Vein©. 

to)   ef*  l~f 


then  m     sin  9 .  n  sin  9. 


dy_    -       .1"!       {j  i-1    £i <2-21) 


dx  m     cos  ©.  n     cos  0 

j-1      *J  i-1      "t 


Equation  (2-21)  could  be  written  in  simple  form, 
dz  -   *Ai   ^ 


sin  9. 
^i 


dx     ri »   cos  ©, 


Z.   cos  g. 

where  A  •  1   if  J    and  9.  are  due  to  poles 

A  -  -1  if  (     and  0.  are  due  to  zeros, 

and    X4  ■  magnitude  of  -8.  +  #. 

©.  *  argument  of  -s. 

Then 


(2-22) 


(2-23) 


In  some  cases  it  is  easier  to  compute  dy/dx  from  the  slope  of  ds/dK 
obtained  by  differentiation  of  -  K  »  p/q  as  shown  in  example  discussed 
at  the  end  of  this  chapter. 
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Fig.  2-3:   Argument  and  magnitude  of  -s   -  6 


K  -  Plane 


S  -  Plane 


1 0  •  IT  +  &<*  ,  ' 


+   4K 


Fig.  2-4:   (a)   Definition  of  K 

t 
(b)   Evaluation  of  S_ 


(b) 
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81 

II -5  -  Graphical  Evaluation  of  S„J 


The  pole  sensitivity  at  any  point  "A"  corresponding  to  a  particular 
value  of  s  can  he  evaluated  graphically  by  summation  of  angles  and  with- 
out evaluation  of  K.  Choose  a  point  C  such  that  AC  is  perpendicular 
to  the  locus  and  evaluate  the  phasor  AC  -  AS   Evaluate  the  value  of 
the  phase  of  K  at  C  by  summing  angles  e  and  let  it  be  ""  +  &<*   This 
procedure  is  based  upon  the  independence  of  the  derivative  of  an 
analytic  function  on  the  direction  of  differentiation.  He  ce  aK  is 
chosen  to  be  A'C  which  is  mapped  into  AC. 
Then:   K  -  j(aot)K 

AS  -   As   -   As   -  S  J  (2-24) 

JK       K 

84 
Therefore  knowing  the  pahsor  As  and  the  angle  A*  one  can  evaluate  S   . 

II-6  -  Curvature,  Pole  Sensitvity  and  Their  Relation 

The  curvature  at  any  point  on  the  root  locus  is  proved  to  be 

(Appendix  II), 

j_  .  .  fa«*>««    .  .    !.<**>„  (2-25) 

z  z 

At  first  glance  it  seems  that  the  curvature  and  sensitivity  are  some!:  ow 
related.  An  example,  discussed  at  the  end  of  this  chapter,  will  reveal, 
however  that  the  curvature  in  that  example  is  constant  while  the  pole 
sensitivity  is  equal  to  cot  <j)  and  consequently  is  a  continuous  variable. 
The  independent  nature  of  the  curvature  and  pole  sensitivity  can  be 
understood  by  noting  Equation  (2-25)  and  the  definition  of  sensitivity. 
Pole  sensitivity  depends  on  first  derivatives,  while  the  curvature 
depends  on  the  second  derivatives,  and  in  general  the  two  derivatives 
are  independent.  As  suggested  by  Ur  we  can  think  of  curvature  as 
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"kinematic"  property  and  the  pole  sensitivity  as  a  "dynamic 

property." 

II-7  -  Open  Loop  Pole  and  Zero  Sensitivity  (Figure  2-5) 

Let  G(s)  be  the  forward  transfer  function  and  L  be  the  loop 
transfer  function.   Then  L  ■  GF  and  the  system  characteristic  equa- 
tion is 

1  +  L  -  0 
and  if  -s.  is  the  system  root,  then 


1  +  L(s) 


8«°8 


where  s  =  s  represents  the  point  on  the  root  locus.   If  K  is  the 
gain  constant  of  L(s)  and  Z  ,  P.  its  open  loop  seros  and  poles,  then 
ve  can  write 


L  -  L(s,K,Zi,Pi) 


and  taking  the  total  differential  of  L: 

n 


m 


dL     -     IL  Jtf     +      \l ,  M    +      2    IiLh?       +      S    \L 


^  dS 


w 


dK 


i«l  6z    dzi 


i-1    bPt  dPi  (2"26> 


On  the  root  locus   1  +  L  ■  0  or  L(s)  »  -1  and  therefore  L(s)  »  constant, 


and  tie  total  differential  dL  is  zero  for  s  ■  -s..     Then  letting 
dL  ■  0  and  8  »   -s,,   Equation  (2-26)   gives, 


18 


8»-8 


d8j  +  "iK 


3 


D"-8 


dK+2 


82, 


8* -8 


^+rf: 


dP, 


8"»-S 


Rearranging  we  get 
1 


ds 


J 


\)B  J   8— S. 


t$)...r;ik):>;K%)A 


B  t  s     itself  is  a  function  of  K,  Z.  and  Ps,  so 
j  i      i 

Sj  -  8j(K,  Z±t   Pj[) 


(2  ?7) 


(2-28) 
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Taking  the  total  differential  of  s.  ve  jet 

K  is,  Jt,     „  lis,  JB      -p    ^s. 


j       bK   K     i  fcZ         i  a?T  V    ' 

Equation  (2-29)  suggests  that  ds.  be  written  as 

fi  s  s 

<»j  "  4J  *   +  I  SZJ  dZl  +  ^  Sp|  dPi  (2'30> 

Equations  (2-29)  and  (2-30)  thus  define  the  sensitivity  of  root  s. 
to  gain  K,  as 

s"J  »  ^  (2-31) 

K 
which  is  the  same  as  Equation  (2-13)  and  the  sensitivity  of  root  s. 


to  open  loop  zero  Z  as 


s'j  £  ±±_  (2-32) 

zi  6zt 

and  the  sensitivity  of  root  s.  to  open  loop  pole  P.  as 

£  £   V  (2-33) 

s .        s  s , 

II-8     -     Relationship  between  S     ,    S  J ,    and  S  J 

Comparing  Equations   (2-27)   and   (2-30)  yields: 

s 


S  * 
K 


s 


K^  H/oK  N  (2-34) 

\  TT'fcS  /  s  -   -s 

W 


i  \Tl  /  s  -  -s 


s 


1  V     bL/iS      /    8    -    -8 


J 


But 


n 


TT  (s  +  Z<) 
L     -     K  J. ^__  (2-37) 

m 

^(8   +    P.) 

1  X 
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Then 


fcK  J  s—  8 


n 

TT  (s  +  z  ) 

1 

j 

1 

(k   )s  -   -8j 


(2-38) 


8    ■     -8 


But  since  1  4-  L(s) 
then  L(s) 


-     0, 


8    «    »8 


J 


1,   and  Equation  (2-38)   yields  to 


8  •   -s 


J 


\W/S    «    -8 


-   I 

K 


Then  (2-34)  becomes 


/  \b*)  s  -   -s^ 


Similarly 


("J  •  -  -. 


n°l 
KIT    (s  +  Z   ) 
1 
m 

"TT  (s+  p.) 

l  x 


s  «  «s 


8   +    Z.       I  -8.    +    Z. 

i    1 8  ■  -s.  j         i 


Then, 


^L/*Z, 
bLAS 


-i/C-Sj  +  zt) 


8    »    «8 


0>L/oS) 


(2-39) 


(2-40) 


(2-41) 


8    ■    -8 


J 


s8j 

KJ 


0R  <?   J      - 

Z 
*i  -s     +  Z1 


(2-42) 


Similarly 


Pi 


K 


s„   -  P, 


(2-43) 


Equations  (2-42)  and  (2-43)  dictate  that  the  sensitivities  to  all 

8  . 

singularities  are  directly  proportional  to  S  J  and  inversely  pro- 
portional to  the  distance  between  -».  and  the  singularity  involved 
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(in  this  case  8  -  -Z.  and  S  ■  -P  are  the  signularities)   Thus  what- 

81  si    8J 

ever  properties  are  found  for  S_J  may  be  extended  to  S_ J   or  S_J . 

A  particular  case  of  Equation  (2-43)  is  for  P  »  0,  i.e.  the 

root  sensitivity  to  the  pole  at  origin  of  the  S-plan. 

8i      s8J 

SpJ  -     SK  (2-44) 

o       (8j) 

81  81 

Then  S  J  is  proportional  to  S   ,  the  constant  of  proportionality  being 

o 

1  ,  a  complex  quantity. 

•j 

II-9  -  Relationship  between  Root  Sensitivity  of  S.  and  Residue 

In  Equation  (2-14),  we  proved  that  in  case  of  unity  feedback, 

SRj  -  -Qj  (2-45) 

where  Q.  is  the  residue  of  T(s)  ■     Kp(s) at  point  -s.. 

J  q(s)  +  Kp(s)  3 

In  this  section  it  will  be  shown  that  if  -s  is  a  single  root,  then 

SKJ  -  F(-8j)  Qj  (2-46) 

where  F(-s.)  is  the  feedback  transfer  function  F  evaluated  at  s  =  -s. 

and  Q  Is  the  residue  of  the  system  transfer  function  at  (-8.).   For 

unity  feedback  and  simple  root  F(-s.)  »  -1  and  Equation  (2-46)  yields 

to  Equation  (2-45).   The  overall  transfer  function  T(s)  of  the  system 

shown  in  Figure  (2-5)  is 

lK8)  1  +  F(s)G(s)       F(s)fl  +  F(s)G(s)]      l*  *'' 

But  L(s)  »  F(s)G(s),  then 

T("  "   F(.)P8|  L(.)1  (2"48) 

'  e  residue  of  T(s)  at  -s.  is 

Q.  -  (.  :  ..)T(.)|        .    (8  !  81)L(S) 

J     |s  -  -8j      F(s)fl  4-  L(s)ls  m     g 
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Then,  let  the  rightmost  expression  he  denoted  as  R.,  i.e.  hy 
definition, 


V-j> 


Q, 


(2-49) 


Then 


R  F(l  +  L)  -  (s  +  s  )L 


Taking  the  derivative  with  respect  to  s  of  both  sides, 

g  [i  +  l]Rj   *n  %  +  m  +  L]  !!i  -   l  +<•+•>  \| 

58 


At  s  -  »s.  i.e.  at  a  point  on  the  root  locus,  L  ■  -1  the  above 

equation  reduces  to 
fcL 


FR 


J 


8  ■  -S 


J 


OR 


v-v  ■  -i/r(-v 

it: 


bs  J  s  ■  -a, 


Comparing  Equations  (2-50)  and  (2°40)  we  get 

8 

s. 


K 


■..)  Qj 


(2-50) 


(2-51) 


Now  turning  to  the  case  of  N  order  root  at  -a.  we  can  expand 


IS 


(.  ♦ ./ 


+  terms  from  other  roots 


where 


rjK 


(N  -  K)! 


+  8 


A 


:i  +  u 


(2-52) 


(2-53) 


Again  defining  the  quantity  inside  the  small  bracket  as  R  ,  then 
repeating  the  operations  as  before,  one  obtains  after  repeated 
differentiation 
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J"      F(-0  \ 

J     N 
8 

8-  -8j 

and  from  Appendix  III 

L  i." J  -  -  —j 

Theorem 

-  N' 

(2-54) 


(2-55) 


When  -s.  is  a  single  order  system  root,  the  sum  of  the 
sensitivities  of  s  to  open  loop  poles  and  zeros  is  equal  to  unity, 
i.e. 

8  8 

2  S7J   +2  SpJ   -  1  (2-56) 

i  Zi      i   Pi 

This  is  easily  seen  by  referring  to  the  construction  of  root  loci. 
If  all  open  loop  zeros  and  poles  are  displaced  by  the  same  amount 
c,    then  all  closed  loop  roots  are  displaced  by  the  same  amount  &, 
that  is  if  dZ.  •  dP  *l    for  all  i,  then  ds  «£  for  all  i.  There- 
fore if  in  Equation  (2-29)  we  put  ds.  »  dZ  ■  dP  «  5  and  dK  »  0  we  get 

h 


A  more  rigorous  proof  is  as  follows.   Rewriting  equations  (2-13) 
and  (2-42)  and  (2-43)  we  get 


i       Zt   -  Bj 

Sf3       -      SK 


i       Sj  -  Pt 
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K 


L  i-l   Zi  "  SJ    1^1  Sj  '  ?i  J 


Then    n   s.     m   s 


Z  spJ   +  2  «  * 


n     S 


in    S. 


■   1 
The  above  relation  is  only  valid  for  simple  poles. 

11-10  -  Graphical  Method  for  Determining  Open  Loop  Pole  and 
Zero  Sensitivities 

Equations  (2-42)  and  (2=43)  suggest  that  to  singularities  which 
are  close  to  system  root  (-s.),  8.  is  more  sensitive  and  for  singular- 
ities which  are  far  away,  8.  is  less  sensitive  until  it  becomes 
insensitive  to  singularities  at  infinity.   This  concerns  the  magni- 
tude of  sensitivities.   But  sensitivities  are  vector  quantities, 
since  the  change  of  a  parameter  may  move  the  roots  in  different 
directions.   It  is  then  helpful  to  make  use  of  Equation  (2-56)  toget- 
her with  Equations  (2-42)  and  (2=43). 

Referring  to  Equations  (2-42)  and  (2-43)  where  Z  and  P  indicate 
open  loop  singularities  and  -8.  is  the  system  root  in  question,  we 

can  draw  a  vector  from  -s  toward  each  and  every  pole,  and  away  from 

J 

each  and  every  zero.   The  length  of  each  vector  is  inversely  propor- 
tional to  the  distance  from  -8.  to  the  singularity  concerned.   Then 
construct  the  sum  U  of  all  these  vectors,  which  is  a  vector  itself, 
(Figure  2-5).   If  U  is  taken  as  a  unity  vector  in  magnitude  and  phase, 
then  the  other  vectors  measure  the  sensitivity  of  s.  to  each 
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S-Plane 


r*-2- 


(a) 


(b)       U 

Fig.  2-5:   (a)  Construction  of  vector  diagram 
(b)  Vector  addition  to  obtain  U. 
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Fig.  2-6:   A  Basic  Feedback  System. 


Fig.  2-7:   Phase  Relationship. 
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singularity  respectively.   So  we  regard  it  as 

U  -  \/o_  -  1  eJo  (2-57) 

and  is  named  by  Thaler  and  Rung  as  "unit-sensitivity"  vector.   This 
unity  scale  does  not  apply  to  sensitivity  to  gain,  S  J.   Actually 
this  technique  is  quite  similar  to  the  graphical  method  of  evaluating 
the  residues  for  simple  singularities,  and  by  putting  U  -  1  /o_   we 
establish  the  necessary  reference  scale.   Phase  of  sensitivity  vector 
is  measured  as  positive  in  the  clockwise  sense  starting  from  the  U 
vector.  This  seemingly  arbitrary  sign  convention  in  fact  comes 
from  Equations  (2-42)  and  (2-43)  which  are  the  basis  of  this  vector 
diagram. 

S«j  -     SK      -       SK        -  SK    (2-58) 


\  .  -  ?1      <*?  -  FT) 


V, 


The  denominator  is  the  vector  from  root  (-s,)  toward  pole  (-P.)  as 
shown  in  Figure  2-7.   The  phase  relationship  of  the  above  equation  is 

8, 


•A 


spJ    *  /  s8J 

(2-59) 


Relation  (2-59)  is  true  no  matter  what  conventions  are  applied  to  the 

measurement  of  the  angles.   Since  the  U  vector  has  been  found  to  be 

81      81  81 

equal  to  1 /o  as  far  as  S  J  and  S."1  are  concerned,  phase  of  S* 

Fi      Zi  Pi 

will  be  measured  starting  from  U  as  zero  phase.  On  the  other  hand 


/    K    and^ ^1 


are  measured  in  the  conventional  way,  i.e. 
starting  from  the  positive  real  axis  and  counting  positively 
counterclockwise. 


J 
In  Equation  (2-59)  the  quantity  /    °K    does  not  depend 

on  i,  i.e.  it  is  the  same  for  all  i's.   Thus  for  each  i, 


Zi 


/,*i  / 

/    P.    ■  constant  -   /   V. 
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/v 

meaning  that  the  larger  the  value  of  /   i   the  smaller  must  be 

/  sSj  /  V 

/   Pi       Since  /    i   is  measured  conventionally  (positive 

/  SSJ 
counterclockwise)  then  /   ?i    must  be  measured  positively  clock- 
wise. 

s  . 
Next  we  prove  that  vector  S  J  lies  on  U.   From  Equation  (2-44) 


s 

Figure  2-7  shows  the  angles 


gives  the  direction  of  U  on  the  S-plane,  it  is  thus  establish- 


s .  s . 

ed  that  on  the  S-plane,  S   always  lies  on  U.   Since  S  J  indicates 

the  direction  in  which  the  root  moves  when  K  varies,  i.e.  the  direc- 
tion of  the  root  locus,  the  above  result  can  be  stated  as:   "At  any 
point  on  the  root  locus,  the  U  vector  is  tangent  to  the  root  locus." 


Also  Equation  (2-59)  reduces  to 


LL--U 


(2-60) 


11-11  -  Locus  of  U  on  the  S-Plane 

Consider  the  plant  to  be  compensated, 

G  ■    .(.  ?\$\  p2)  <2-6l> 

and  the  desired  dominant-root  location  -s.  as  indicated  in  Figure 

J 

2-8.  A  spirule  measurement  shows  that  an  additional  phase  of  +  (p 

is  needed  at  location  -s . .   Thus  a  lead  network  with  a  zero  at  Z  and 

J 

a  pole  at  P  is  needed.  The  question  is,  how  does  the  tip  of  the  U 


33 


vector  move  on  the  S-plane,  when  Z  and  P  take  all  possible  values 

on  the  negative  real  axis? 

Now  sensitivity  vectors  should  be  drawn  for  all  the  plants  poles 

and  zeros,  ignoring  for  the  moment,  the  need  for  compensation.  All 

these  vectors  can  then  be  added  to  obtain  the  vector  sum  QI,  which 

may  be  called  the  "uncompensated"  unity  vector.  Then  the  geometric 

locus  of  U  is  a  circle  with  a  center  I  and  a  radius  of  r  =    1     . 

d  sin  i 
This  circle  is  named  the  (U)  circle. 

Proof; 

QI  would  be  the  unity  vector  without  "P"  and  "Z"  (Figure  2-8a). 

A        . 

When  P  and  Z  are  added  so  that  PQZ  =  <p,  a  sensitivity  vector  must 

be  assigned  to  each  of  these  singularities,  Qm  =  1  is  assigned  to 

QP 
the  pole  and  Qu  =  1  is  assigned  to  the  zero.   The  vector  QI  then 

QZ 
increases  by  the  vector  quantities  IM  =  Qm  and  MU  =  Qu.   QU  is  the 

final  U  vector.  Note  that,  in  Figure  2-8a,  IN  and  MU  are  two  equal 

but  opposite  vectors. 

Figure  2-8b  shows  that  when  P  moves  along  the  real  axis,  since 

Qm  =  _1  ,  m  moves  on  a  radius  of  r  =  _1,  because  of  inversion  of 

QP  2d 

real  axis  with  Q  as  center  and  ratio  of  inversion  one.   Such  a  circle 

is  referred  to  as  the  (m)  circle.   If  we  translate  this  (m)  circle 

by  QI,  then  we  get  (M)  circle  which  is  the  locus  of  point  M.   So  the 

locus  of  point  M  is  the  (M)  circle,  or  radius  R  =  _1  and  with  I  as 

2d 
its  upper  most  point.   Figure  2-8c  shows  that: 

IU  =  MN  -  2R  sin  MON 

2 

IU  =  2R  sin  <|>  =  I   sin  §   ■  constant      (2-62) 
d 
Hence  U  moves  on  a  circle  centered  at  I  and  of  radius  r  =  1_  sin  <^. 

d 
So  the  following  result  can  be  stated: 
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U  circle 
radius  (1/d)  sin 


(M)  circle 
radius  1/d 


(m)  circle 
R  1/d 


Fig.  2-8:   Geometrical  Locus  of 
Point  U. 


N 
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"The  geometric  locus  of  the  tip  of  the  unity  vector  is  a  circle, 

centered  at  I  and  of  radius  r  =  1  sin  If,  where  I  is  the  tip  of  the 

d 
'uncompensated'  unity  vector,  d  the  imaginary  part  of  the  dominant 

system  root  and  (J)  the  phase  shift  to  be  introduced  by  the  compensation." 

Note  that  there  are  limits  to  the  geometric  locus  of  U  on  the 

(U)  circle.  One  limit  corresponds  to  the  extreme  case  where  P  »  «>, 

the  other  is  for  Z  ■  0.  Later  it  will  be  shown  how  to  obtain  these 

limit  points  on  the  (U)  circle. 

11-12  -  "Use  of  the  U  Locus" 

Starting  with  the  uncompensated  plant  and  a  desired  location  for 
the  dominant  root,  it  is  always  possible  to  be  able  to  draw  the  (U) 
circle  and  be  able  to  obtain  the  two  limit  points  on  the  U  locus  on 
the  circle. 

Conversely  after  a  point  U  has  been  chosen  on  the  U  locus,  it 
is  possible  to  perform  the  graphical  construction  in  reverse  order 
and  thus  drive  the  corresponding  Z  and  P,  i.e.,  obtain  the  compensation 
needed.   The  following  will  explain  how  to  choose  the  U  point  on  the 
U  locus,  which  is  the  seence  of  this  construction  technique. 

In  Figure  28a,  note  that  OJ  is  perpendicular  to  IU  where  J  is 
the  midpoint  of  arc  MN.   P  and  Z  can  be  obtained  by  drawing  OJ  per- 
pendicular to  IU,  cutting  the  (M)  circle  at  J,  drawing  arcs  JM  ■  JN  ■ 
<J>.   P  is  then  determined  by  drawing  OP  parallel  to  IM,  and  Z  is  det- 
ermined by  drawing  QZ  parallel  to  IN. 

An  example  is  given  in  Figure  2-9a,  where  a  desired  U  point  is 
shown  (an  arbitrary  choice  for  illustrative  purposes).  OJ  is  drawn 
perpendicular  to  IU.  Next,  arc  JM  ■  arc  JN  »  <j)  is  drawn.   In  Figure 
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(a) 


*     U  circle 


Z 

(b) 

/"N 

t     I    i 

Z 

it 

'v 

p         ly 

-pi    ' 

j* 

Fig.  2-9: 

Determination 

of 

P  and  Z  for 

given  U. 
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(•) 


(b) 


U  circle 


Fig.  2-10:   Determination  of  the  limit  point*  of  geometric 
locus  of  U. 
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2-9b  P  is  obtained  by  drawing  QP  parallel  to  IM  and  Z  is  obtained  by 
drawing  QZ  parallel  to  IN. 

Now  the  problem  recurs  of  determining  the  limit  points  of  the 
U  locus,  using  the  same  construction  procedure.   The  extreme  right 
limit  point  corresponds  to  the  case  where  one  compensator's  singu- 
larity (here,  the  zero)  is  at  the  origin.   The  extreme  left  limit 
point  corresponds  to  the  case  where  the  other  singularity  is  at  -  oo„ 
Figure  2 -10a  shows  how  to  obtain  the  right  limit  point,  U    Draw- 
ing IN  parallel  to  QO,  because  Z  is  now  at  0  and  take,  arc  N  Js  f 
Drawing  IU  perpendicular  to  OJ  ,  then  results  in  U  . 

Figure  2-10b  shows  how  to  obtain  the  left  limit  point,  U. . 

Draw  IM.  parallel  to  QP,  which  is  horizontal  since  P  is  at  -  oo,  and 

take  arc  M..J.  =  <j).   Drawing  IlL  perpendicular  to  OJL  ,  then  results 

in  U1 .   The  geometric  locus  of  U  is  the  part  of  (U)  circle  between 

IL  and  U  . 
1      r 

11-13   -  Design  Techniques 

In  the  preceeding  section  it  has  been  derived  a  method  for 

finding  P  and  Z,  given  the  desired  location  of  point  U  on  its  locus, 

i.e.  the  magnitude  and  phase  of  QU,  the  unity  vector.   This  method 

is  restated  below  in  a  step  by  step  form,  Figure  2-8a. 

Step  1  -  Considering  Q  as  if  it  already  were  a  point  on  the  root 

locus,  draw  the  vector  diagram  at  Q  for  the  uncompensated 

system  and  obtain  QI,  the  "uncompensated"  unity  vector. 

Step  2  -  Draw  the  circle  diagram,  composed  of  the  (M)  circle,  of 

radius  R  =  _1  and  whose  uppermost  point  is  I,  and  the  (U) 

2d 
circle  centered  at  I  and  of  radius  r  ■  1^  sin  <j>.   Fix  the 

d 
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limits  of  the  locus  of  U  on  the  (U)  circle. 

Step  3  -  Given  QU  as  the  desired  unity  vector  for  the  compensated 
system,  draw  OJ  penperpendicular  to  IU  which  cuts  (M) 
circle  at  J.  Draw  angles  JON  =  (|>.  Then  Z  and  P  are  deter- 
mined by  drawing  QZ  parallel  to  IN  and  QP  parallel  to  IM. 
In  this  section,  it  will  be  seen  how  the  "U"  vector  (i.e.  the 

location  of  U  on  its  locus)  is  selected  to  satisfy  a  particular 

condition. 

(a)  Design  for  Minimum  Root-Sensitivity. 

Since  the  U-vector  is  the  scale  used  to  measure  the  individual 
sensitivity  vectors,  the  larger  the  scale,  the  smaller  the  magnitude 
of  the  sensitivity  measures.  Thus  one  possibility  is  to  design  the 
compensation  for  maximum  magnitude  of  the  U  vector,  that  is  minimum 
root  sensitivity  to  open  loop  singularities.  Figure  2-8a  shows 
that  maximum  magnitude  of  QU  is  obtained  by  placing  point  U  near  the 
lower  most  part  of  the  (U)  circle  or  more  exactly  on  the  extension 
of  QI.   If  such  point  is  not  within  the  locus  of  U,  then  it  can  not 
be  a  location  for  U  and  one  must  select  the  lowest  point  which  is  on 
the  locus.   In  Figure  2-9a,  this  is  point  U  .  Thus  QU  is  the  selected 
unity  vector,  and  with  this  given,  one  can  proceed  to  the  3  step  pro- 
cedure outlined  in  the  beginning  of  this  section.  With  such  a  design 

s  .      s  . 
S  J  and  S  J   are  all  minimum,  for  all  i's. 
i      Zi 

(b)  Design  for  Constant  Damping  When  K  Varies. 

Another  practical  problem  is  to  compensate  a  system  in  such  a 
way  that  when  gain  K  varies  about  its  nominal  value,  the  dynamic 
response  of  the  system  does  not  change.  This  calls  for  a  constant  3, 
i.e.,  a  root  locus  that  remains  tangent  to  the  radial  line  OQ  at  the 
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neighborhood  of  Q,  Figure  2-8a.   How  can  Z  and  P  be  found  to  obtain 

such  a  root  locus?   It  is  now  shown  that  this  can  be  done  by  merely 

selecting  point  U  so  that  the  unity  vector  QU  goes  through  the  origin, 

0,  of  the  S-plane.   In  other  words  choose  U  so  that  Q,  0  and  U  are 

in  line.   If  the  locus  of  U  does  not  permit  such  a  choice,  this  means 

it  is  not  possible  to  obtain  a  constant  T  about  Q  for  the  given  system. 

One  can  then  choose  the  best  solution  available,  by  taking  the  "U" 

location  that  is  closest  to  a  straight  line  with  Q0. 

The  above  statement  can  be  proven  very  simply  if  one  recalls 

Equation  (2-44)  that  s. 

s.        S  J 

SpJ   -  —p-  (2-63) 

o        j 

s  . 

The  specification  here  is  to  force  S  J  to  have  the  same  direction  as 

K 

0Q.   Thus  making  s.  move  on  a  radial  line  when  K  varies.   This  means 


L 


J 


S  J    must  be  equal  to  the  phase  of  Q0  which  is  also  phase  of  s . . 


But  from  Equation  (2-63) 

-  0         (2-64) 


La.-  La.  ■  l 


ii 


The  phase  of   j  must  be  zero,  this  means  that  the  sensitivity  vector 

s .  ro 

Sp-*  must  be  on  the  U-vector  or  conversely,  the  U  vector  must  pass 

through  P  at  the  origin  of  the  S-plane. 

s  . 
A  faster  way  to  prove  this  is  that  we  show  that  S  J  always  lies 

e  si 

on  the  U  vector.   In  order  to  keep  3  constant,  S   must  be  radial, 
thus  U  vector  must  be  radial. 


11-14  -  Limiting  Behavior  and  Special  Cases. 

The  magnitudes  of  gain  sensitivities  can  cover  the  entire  range 
of  values  from  minus  to  plus  infinity.   Yet,  intuitive  notion  of 
"sensitivity"  as  a  general  concept  in  closed  loop  systems  makes  part 
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of  this  range  unreasonable.  One  part  of  the  problem  is  a  direct 

consequence  of  the  sensitivity  definition,  while  another  is  associated 

with  its  first  order  approximation  nature.   A  better  understanding  of 

both  facets  can  be  gained  by  an  examination  of  limiting  cases. 

In  general,  closed-loop  poles  depart  from  open-loop  poles  for 

low  values  of  gain,  and  proceed  to  either  open  loop  zeros  or  unbounded 

values  as  the  open  loop  gain  becomes  very  large.  The  gain  sensitivity 

as  given  by  Equation  (2-13)  could  be  written  in  the  form 
s  . 

SRJ   =  1 (2-65) 

n      ,       m+n     . 

77  1 +  yr   1 

f-*.   Z,  -  s.    f\      s.  -  P. 
i=l   i    j    i=l   j    i 

As  K  approaches  zero,  the  closed  loop  root  (-s.)  approaches  the  open 

loop  pole  from  which  it  derives,  i.e.  s.— >P. .  Then  the  term,  the 

above  relation  reduces  to 
s 


S  J 
SK 


K-»0 


UVJ 


J+     0      (2-66) 


Similarly  as  K  becomes  very  large,  n  of  the  closed-loop  poles  approach 

open  loop  zeros.   If  the  jth  closed-loop  pole  is  one  of  these,  and 

it  approaches  the  ith  open  loop  zero  so  that 

s. >  Z. 

J       i 

Then,        s .  ] 

sv3  »      1  >  0      (2-67) 

*  J  K-»co       /    1      \ 

C  zi  -  sj  ; 

Finally,  m  of  the  closed-loop  poles  have  no  zeros  to  go  to,  and  hence 


become  very  large  relative  to  P  and  Z..   The  sensitivity  for  these 


poles  is 


s  . 
J 


S  J      >  m  _l_  =  _1L    (2-68 

s.  ^  Z.,P.        TT    8,         m 


i-1 
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Here  m  =  (order  of  denominator)  -  (order  of  numerator).  When  the  gain 
is  sufficiently  large  for  the  open-loop  zero  db  line  to  intersect 
the  high  frequency  asymptote  on  Bode  plot,  the  open  loop  transfer 
function  is  approximately 

L(s)  = 


sm 

where  m  =  (order  of  denominator)  -  (order  of  numerator).   Then  for 
the  points  on  root  locus 

1  +  L(s)|        =  0 

Is  ■  -s , 

1 

or  1  +     K      =  0 

(-Sj) 
or  s.  -   -   V^K"  (2-69) 

Thus,  the  sensitivity  of  the  unbounded  pole  from  (2-68)  will  be 


•'1    — 

jsJ>>Z1,P1 


V  I      *      ~       "^  -B 

m 


Equation  (2-70)  indicates  that  sensitivity  increases  as  the  m  root 
of  K,  as  gain  is  increased,  although  for  the  finite  gains  the  sens- 
itivity is  always  finite.   Referring  to  Figure  2-11,  the  above  special 


cases  can  be  simplfied  as:  /*  s^ 


For  poles  going  towards  "O.L.  zeros"  <  „    i  ^  Q 


S„J  -♦  0 


For  poles  going  towards  infinity      [   S  J  =  _j_       in  general 

m 

svl      ~      ~   SZ  ~k  wnen  L  =  K 

m  m 

s 

Another  circumstance  in  which  the  sensitivity  can  become  very 

large  is  revealed  by  Equation  (2-45),  and  that  happens  when  Q.  (the 
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Fig,  2-11:   Poles  going  towards  open 
loop  zeros. 


-8 


•P    -z 

ri    i 


Fig.  2-12:   Poles  going  towards  infinity. 
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residue  evaluated  at  pole  -s  )  Is  not  finite.  This  is  to  be  expected 
since  the  sensitivity  factors  defined  thus  far  have  not  considered 
multiple-order,  closed-loop  roots.  As  far  as  the  gain  is  finite, 
an  infinite  gain  sensitivity  always  indicates  multiple-order  closed- 
loop  poles. 

A  special  situation  of  considerable  interest  can  occur  when  a 
closed  loop  root  lies  between  an  open-loop  pole  and  zero  which  are 
much  closer  to  each  other  than  to  all  other  open- loop  poles  and  zeros, 
This  is  the  so  called  dipole  case  as  in  Figure  2-12.  The  sensiti- 
vity for  the  bounded  closed-loop  pole  will  be,  approximately 


v  -  . 

1 

1 
8j-Pi 

-  <Vi><VFi> 

1       . 

zi-pi 

Z     -  s 
i         j 

(2-71) 
Differentiating  equation  (2-71)  to  find  values  of  s,  which  make, 

K 


81 
sensitivity  maximum,  reveals  that  maximum  value  of  S„J  will  occur 


Z  +  P  81 

when  s  ■  i    1  for  which  S„  becomes 


SK 


-  1  (Z.  -  P.)  (2-72) 

max 


11-15  -  Sensitivity  at  Irregular  Points 

Expanding  the  gain  K  ■  K(s)  in  series  to  include  the  higher 
order  terms  we  have 

K  «  K  +  dK  As  +  1  d*J|  As2  +  -— -  (2-73) 

°   ds     2.2 
ds 

and  then  assume  the  first  n-1  derivatives  vanish.  This  means  then, 

AK-K-K  -  1   d^As" 
o      I    " 
n.   .  n 
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or 

As     -  n  /n!  Ak_      -    n  /   n! 

V^K     (2-74) 

Since  As  is  finite  it  is  suggested  to  be  substituted  for  S_J 

at  irregular  points. 

For  any  AK  >  0,  ds  has  n  values  separated  by  2TT  which  represents 

n 
the  outgoing  branches  in  Figure  2-13.   If AK  <  0,  n  branches  with  the 

same  separation  as  before  are  again  obtained,  but  they  are  displaced 
from  the  previous  set  by  TT/n,  that  is  the  two  set  intersect.  This 
could  be  stated  as:  The  incoming  branches,  which  at  their  junction 
represent  an  nth  order  closed-loop  pole,  are  evenly  spaced  and  sep- 
arated from  each  other  by  2lT/n.   The  outgoing  branches  are  also 
separated  from  each  other  by  277/n  and  are  midway  between  the  incoming 
branches . 

11-16  -  Sensitivity  Functions  for  Alternate  Transfer  Function  Forms 

The  transfer  function  could  be  written  either  in  root  locus  or 

Bode  form.  For  Figure  2-6,  open-loop  transfer  function  can  be  written 

as  n  .    n-1  , 

s  +  a,s    + +  a 

n 


c<.>  -  gjsi  -  k  -zr± 

x  ■     E(s)        trt-m 


m+*n  - 1  .  . 

8     +  bi§  +  ....  bmfn 

2  a.."'1  S  (•  +  Z.) 

K  i-0  m   R  i-1 L      (2-75a) 

mfn     ,   .  mfn 

2  b.s^"-1  TT  (s  +  p  > 

i-0  x  i-1 


m  / 

:•  i-i  vzi 


+  i 


mfn  /  s_   -> 


(2-75b) 
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The  Equations  (2-75a)  and  (2=75b)  represent  the  root  locus  and  Bode 
form  respectively. 

The  equations  for  the  gain  sensitivity  are  the  same  whether  open- 
loop  transfer  function  is  in  root  locus  or  Bode  form,  i.e. 


"j 


K 


s 


(2-76) 


This  is  because  in  comparing  (2-75a)  and  (2~75b)  it  is  obvious  that 

„.       TT 


K  i  i 


K  x  constant 


i  i 


and  then, 


bs 
K 


*s 


K' 


It  is  however  necessary  to  modify  the  open  loop  pole  and  zero  sensi- 
tivities for  terms  which  are  written  in  Bode  form.   For  poles  and 
zeros  which  appear  in  open  loop  transfer  function  in  Bode  form,  i.e. 


&> 


or 


§:- 


,  the  sensitivities  are 

s 


s 


j 


3,  S, 


j 


zi(zi  -  Sj) 


and 


P. 


(2-77) 


(2-78) 


Frequently  open-loop  zeros  and  poles  will  occur  as  complex 
conjugate  pairs,  and  variations  in  the  system  will  change  both  zeros 
or  poles.   For  example,  consider  a  complex  pair  of  zeros,  Z  and  Z«, 
which  are  defined  by  their  frequency,  (so,   and  damping  ratio  ^>  *••*• 
Z,  «  Sfc)  +  jG>Vl  -  ^2 


1 

z2  -  SP  -  j"/i  -S2 


(2-79) 
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For  this  situation  it  is  useful  to  define  frequency  and  damping 
ratio  sensitivities  as 

<j  -  ii   -  i!i    fi   +  lit    fa 

5         V£         bzx      d§  Tz£      d^s 

Where  ^  and  £  are  given  according  to  (2-42)  or  (2-43)  (In  case  of 

azi     *Z2      dZx   dZ1  dZ2      dZ 
polee).  Components  as  ~"T£j""»  Tei  7£T  »  *nd  tt  can  easily  be  eval- 
uated from  (2-79).  Then  when  the  complex  pair  of  zeros  appear  in  open- 

2      2 
loop  transfer  functions  in  root  locus  form,  i.e.  (s  +2^<i*+K>),  the  fre- 


«.<o.  are    .      am.   ,  S  ••  .  ..   ■  r _  (2-81) 


i***t] 


8_  , 

10 


quency  and  damping  ratio  sensitivities  for  a  complex  pair  of  poles  or 

Sj  -  236>B   +6*  8*  -  2  $6)8   + 

Where  the  upper  sign  is  to  be  used  for  zeros  and  the  lower  one  for  poles. 

If  the  term  is  written  in  Bode  form,    i.e. 

sitivities  are         "j      £fl  (ft°    "  81)SK^ 

^°     '     W  s2-25*sj+W2 

ij  -    +  !!f  SK3  (2-82) 

s2  -  2^8j  +  *2 

For  some  cases  it  may  be  more  convenient  to  define  a  complex  pair  in 
terms  of  their  real  and  imaginary  parts,  i.e. 

Zx     -  a  +  jb        Z2  •  a  -  Jb 

2  2    2 

In  this  case  with  the  terms  in  root  locus  form  [s  +  2as  +  (a  +  b  )], 

the  sensitivities  are 

s,     *>s      Is  dZ.     ^s      dZ 

c  J   m       —J,  m      __J.  m      __£   4.   1  m      __£ 

a      ^a     ^Z   p   da  ^  ^Z2  •   da 

8  8 

s  s 

-  £f  '       1 

I  (a  +  jb  -  Sj)   +   (a  -  jb  -  8j) 


48 


2(a  -  Sj)  SKJ 


2      n  2      ,  i 

s.   -  2as.  +  a  +  b 


and  in  a  more  general  form  as  in  case  of  Equation  (2-81)  it  could  be 
written  as  s, 


(2-83) 


8j  .  i!i 

a             ~ba 

8j  .  ^i 

m 

as 

+2(a  -  8j)  SRJ 

s.      -  2as.  +  a    +  b 
+  2b  S® j 

b               ib 

2                        2         2 

8j     '-"j 


where  the  upper  sign  is  to  be  used  for  seres  and  the  lower  one 

/  2   2    2 
for  poles.   If  the  term  is  in  Bode  form  f  [s  /(a  +  b  )]  + 

2    2      \ 
[2as/(a  +  b  )]  +  lj  ,  the  sensitivities 


are 


';  -  f  t  », 

i  2         2 

a*  +  bz 


±2V 
a2+b2 


(2=84) 


11-17  -  Example 

In  Figure  2-14,  suppose  the  gain  ¥.,    is  given  as 


2 

s  +  1 

8 


1 
>8  -  — 

S 


Find  radius  of  curvature,  and  pole  sensitivity. 
We  can  write 


dK 
ds 


2 
dTK 

a.2 
A 

.  3 

ds 


6 


Double  points  occur  at 


or  S  «*  +  1.   Substituting  S 


1«*\ 


ds 


Fig.   2-13:     Outgoing  Branches   (n  -  2) 


S  -  X  +  jY 


Fig.    2-14:      Example  of 
Section  II 
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we  get 


dK 

■ 

=•1   + 

1 

ds 

e2j0 

■ 

=  1   + 

.-2J* 

■ 

2Ja" 

**  8j*  i 

e^ 

2j 

m 

-2je 

"^  sin  $ 

s 

2  sin 

^e^  + 

2 

On  circular  portion  of  the  root  locus, 
d  d~ 

Rcd7  2 

The  pole  sensitivity  is 

8i       ^  2  J  ,       2 

o  J       w  ds  S   4-  1        8 


K       dK  s    w  ,    2 

1  -  8 


1 

8  -  — 
8 


On  the  circular  portion 


81      Si     ej*  +  e°j'         \i  j 

V  "  •       ji      °j4    "   °je    cot ' 

The  pole  sensitivity  at  a  ■  -1  (a  double  point)  is  infinite  and  hence 
the  suggested  substitute,  i.e.     As   is  evaluated.   Here  n  ■  2  and 

As 


Yak 


a»    4-  1  evaluated  at  8  ■  -1 


2! 
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Ill 

SENSITIVITY  AMD  LARGE  PARAMETER  VARIATIONS 
III-l   -  Discussion 

This  section  deals  with  the  most  practical  aspect  of  the  sen- 
sitivity and  is  quite  utilized  in  design  consideration  of  the  feed- 
back control  systems.  It  discusses  sensitity  in  two  particular  cases 
one  with  leakage  transmission,  and  the  other  case  when  leakage  trans- 
mission is  zero.   It  relates  the  sensitivity  of  the  system  to  return 
difference  and  null  return  difference. 

In  previous  sections  it  has  been  discussed  that  in  moat  cases 
we  are  facing  the  variation  in  the  values  of  an  element  in  an 
engineering  system.  This  element  could  be  a  vacuum  tube,  transistor 
or  sometimes  could  just  be  a  passive  element.  It  la  possible  by 
means  of  feedback  around  the  troublesome  element,  to  achieve  the 
desired  reduction  in  sensitivity.  Let  the  troublesome  parameter  be 
k,  as  in  Appendix  1.  We  relate  k  to  the  controlled  source  S  and 

control  variable  C  i.e.  S  »  kC  and, 

kt  ,t 

T  -  t   +  cl  os  (3-1) 

1      oi  +   1  -  kt  lJ  l) 

C8 

with  the  fundamental  signal  flow  graph  as  shown  in  Appendix  1. 
We  defined  the  classical  or  system  sensitivity  as, 

T  m     AM_  .   >T/T  (3.2) 

T 
so  S,  is  the  relative  change  in  T  divided  by  the  relative  change  in 

k  for  infinitesimally  small  changes  in  T  and  k  only. 

Case  1; 

Consider  the  case  v/here  t  .«  T,  i.e.  leakage  transmission  is 

negligible.   This  is  a  situation  in  many  feedback  systems.   By  straight 
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forward  differentiation  of  Equation  (3-1)  we  get 


and  from  that 
Sk  a 


WT 

Tk/k 


k 


kt  .t 
ci  OS 

kt 


cs 


k<x  ;  ktcs> 

kt  ,t 

Ci  OS 


kt 


cs 


k  \T 

T  Tk~ 


t  .t  (1  =  kt  )  +  t  (kt  ^t  ) 

ci  os       cs     cs   ci  os 

(1  -  kt  )' 
x      cs' 


t  „t 

ci  OS 


[l  -  kt  ) 
cs 


for  t 


oi 


(3-3) 


where  F,  is  the  return  difference  for  reference  k.   It  is  clear  now 
k 

why  it  is  insisted  that  all  t  „'s  in  the  fundamental  feedback  equation 


should  be  independent  of  k.   To  have  a  small  sensitivity  then 


F, 


should  be  large,  i,e.  (°kt  )  which  is  usually  called  the  loop  gain 

should  be  a  large  number. 

T 
Since  t   and  possibly  k  is  a  function  of  frequency,  so  is  S.  . 

CS  K 

T 
and  one  must  consider  the  frequency  range  over  which  a  small  S.  is 

desired.   It  is  interesting  to  note  that  as  the  magnitude  of  the  loop 

transmission  determines  the  sensitivity,  its  phase  also  determines 

the  stability  of  the  system. 

Case  2: 

When  the  leakage  transmission  t  »  4   0,  then  by  straight  forward 

ol 


differentiation  as  in  Case  1  we 
&T 


t  .t 
ci  os 


TT 


(1  -  kt   ) 
v      cs 


and  S 


k 


k  *T 
T  Tk 


kt  „t 

ci  os 


T(l  -  kt 


cs 


S3 


Pig.  3-1:   Signal  Flow  Graph  for  Example  (1). 


Input 


1 


Fig.  3-2:   Signal  Flow  Graph  for  Example  (2). 


1 

»■  9 


Input  1      a 


1     output 


Fig.  3-3:   Signal  Flow  Graph  for  Example  (3) 
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This  could  be  written  in  the  form 

1 


T 

sk  - 


kt 


cs 


and  since  from  Equation  (3°1 
T  « 


kt  ,t 
ci  OS 


oi 


kt 


cs 


then 


1 

1  -  kt 


kt  „t 
ci  os 


:i  -  ktcs> 


F 


k 


(3-4) 


S,   can  also  be  written  in  terms  of  return  difference  and  null  return 
k 

difference.   From  Appendix  1,  Equation  (A3»3),  we  have 


oi 
T 


F' 


then 


,T 


K 


(3-5) 


where  F,  is  the  return  difference  with  reference  to  k  and  F'  is  the 
null  return  difference  also  with  reference  to  k.   Equation  (3-4) 

suggests  a  practical  mean  to  measure  the  sensitivity.   It  is  noted 
that  in  case  t   «  0,  Equation  (3-4)  reduces  to  Equation  (3-2).   The 
applications  of  Equation  (3°4)  are  illustrated  by  a  number  of  simple 
examples . 
Example  °  1 

In  the  single°loop  system  Figure  3=1,  if  the  sensitivity  desired 
is  with  respect  to  the  forward  transmission,  the  leakage  transmission 
is  zero,  because  with  K  ■  0  there  is  no  transmission  from  input  to 
output.   The  sensitivity  in  this  case  is  the  reciprocal  of  the  return 
difference,  i.e. 


Example  -  2 

In  the  same  single-loop  system,  If  the  sensitivity  of  interest 
is  with  respect  to  the  feedback  transmission, ft  ,  the  leakage  trans- 
mission will  then  be  K,  and  Equation  (3-5)  yields 

T        1     /        K      \      0IL 
Sf  "   1  -M  \  "    K/(l  -^K)  /  "   1  -  £K      (3-7) 

In  this  example  it  is  noted  that  when  0-  0,  i.e.,  there  is  no  feed- 
back around  k,  then  the  leakage  transmission  t.»lxKxl»K. 

oi 

Comparison  of  Equations  (3-6)  and  (3-7)  demonstrates  the  well  known 
fact  that  a  value  of  K  large  with  respect  to  unity  results  insensitive 
to  changes  in  K,  but  with  the  closed-loop  gain  essentially  equal  to 

1  (when  K  is  very  large)  the  sensitivity  of  T  with  respect  to  I 

'  *     o 

approximates  the  value  of  -1  ■  1  /18  . 

Example  -  3 

Figures  3-2  and  3-3  represent  alternative  designs  for  a  feed- 
back amplifier.  Three  amplifier  stages  are  available,  each  with  a 

gain  a,  and  the  overall  gain  to  be  realised  is  specified  (less  than 

3 
a  ) .  Negative  feedback  is  to  be  introduced  to  reduce  the  sensitivity 

of  the  overall  gain  to  variations  in  supply  voltage.  To  a  first 

approximation  this  voltage  variation  can  be  considered  equivalent  to 

a  variation  in  each  value  of  a.  The  question  to  be  answered  is  which 

configuration  gives  a  lower  value  of  sensitivity  of  the  overall  gain 

with  respect  to  a?  For  Figure  3-2,  the  overall  gain  T.  is 

,3 

T, 


H  "    [t*t] 


Th 


e  sensitivity  of  T  with  respect  to  a  is  three  times  the  sensitivity 
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of  an  individual  stage  with  feedback,  or 

s.1  ■    i  -v  <3-8> 

For  the  second  configuration,  with  a  single  overall  feedback  path, 

the  overall  gain  T-  is  3 

2         „         a 

T    s  ^ 

2    i  ■  ?/ 

and     ST2  .  3    .  (3.9) 

1  -  v 

The  two  sensitivites  (3=8)  and  (3-9)  air©  to  be  compared  on  the  basis 

2 

of  the  equality  of  T.  and  T  .   This  equality  establishes  the  relation 

(1  -  /^a)3  .  1  -  ^a3  (3-10) 

substituting  (3-10)  into  (3-9)  we  get 

T 
„  2  3 


ci  -  far 

Thus,  the  single  overall  feedback  path  results  in  a  system  with  less 
sensitivity,  to  changes  in  a. 

Ill -2  -  Large  Parameter  Changes. 

A  serious  shortcoming  of  Equations  (3=1),  (3=2),  and  (3-3)  is 

that  they  apply  only  for  inflnltisimally  small  changes  of  k.   There 

is  therefore  uncertainty  as  to  their  applicability  for  moderate  or 

large  changes  in  k.   It  is  therefore  found  useful  to  use  a  new 

definition,  defined  as  follows: 

Let  T  ,  k  represent  the  nominal  or  original  design  value  of 
o   o 

the  system  transfer  function  and  of  the  element  under  consideration 
respectively,  and  let  T  ,  k  be  the  corresponding  (final)  values 
at  the  new  value  of  k„   Thus, 
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I£  . 

T    +  AT 

kf  ■ 

k    +  Ak 
o 

(3-11) 

The  new  sensitivity  function 

is  defined 

as 

T     D 
Sk     * 

AT/T 
Ak/kf 

(3-12) 

Case  1: 

When  it  is 

assumed  that 

'oi     -      °« 

then  from  Equation  (3-1), 

T       - 
o 

kt   .t 

O    Ci    08 

(3-13) 

1  -  k  t 
o  cs 

if  - 

krt   .t 
£  ci  os 

Tf"To 

(3-14) 

1  -  k^t 

£  cs 

AT/T 
Ak/kf 

Sk     - 

Tf 

(3-15) 

kf-ko 

kf 

Substituting  (3-13)  and  (3-14)   into  (3-15)  we  get 

l-To 
Tf 

1 

k  (1  -  k-t     ) 

o             f  cs' 

Sk    " 

k,(l  -  k  t     ) 
fv           o  cs7 

k     -  k 

f         o 

k£"ko 

kf 

« 

1 

1 

1 

(3-16) 

1  -  k  t 
o  cs 

Fk 

1  +  L 

0 

where  F.   is  the  initial  return  difference  and  L  is  also  the  initial 

ko  ° 

loop  transmission.   It  is  thus  shown  that  when  the  leakage  transmis- 
sion t   is  zero,  the  classical  and  new  definition  have  the  same 
oi 

value,  although  they  are  defined  differently.  With  the  new  definition 

there  is  no  uncertainty  as  to  the  effect  of  large  parameter  variation. 

Equation  (3-16)  is  exact  only  if  t   is  zero,  but  it  may  be  used 

oi 

whenever  t  «  Tf ,  T  .  Fortunately  t  .  is  zero  for  a  very  large  class 
of  feedback  amplifiers.  From  Equations  (3-11)  and  (3-16)  we  have 
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T     -  T  t 

f          o  ,          o            AT            ^T     Ak 

s  »      1    «  — —  ■     a —     *     S 

T  T              T                 k     k 


k  t 

O    Ci 


k  t 
o 


Therefore 


_     »     1 


a  -  j*  >  (  j-rh-  ) 

f  \  O    C8      ' 

\ri  ko"cs  )\    1-  kotcg    / 


k 
and  _o 

T  k 


+  L 


t!   '       "TTT—  <3"17> 

f  o 

where  L  ■  -k  t   is  the  original  or  the  nominal  loop  transmission 
o     o  cs  r 

and  Equation  (3-17)  relates  the  original  and  final  values  of  the 
closed  loop  transfer  function  which  is  of  very  practical  design 
application. 
Case  2 

In  previous  cases  we  assumed  t  .  was  negligibly  small  in  com- 

oi 

parison  with  T  and  T  over  the  range  of  variation  of  k„   This  is  not 

always  true.   Suppose  in  a  system  whose  signal  flow  graph  represent- 

ation  is  that  shown  in  Figure  3-4  parameter  D  varies  substantially, 

and  feedback  is  to  be  used  to  reduce  the  system  sensitivity  to  the 

variations  in  D.   If  we  choose  D  e  k  of  the  fundamental  feedback 

equation,  then  t   may  not  be  small.  However  the  problem  can  be 

handled  by  method  in  Case  1,  if  we  simply  let  k  represent  the  entire 

element,  i.e,  fc  u       A(B  +   ED) 

1  ■ 


Fig.  3-4:   A  Signal  Flow  Graph  with  Leakage  Transmission. 


R                      E 

G(s) 

C 

J 

i 

H(s) 

Fig.  3-5:   Basic  Feedback  System. 


G(s) 


Fig.  3-6:   Basic  Unity  Feedback  System. 
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The  element  k  is  then  incorporated  into  the  overall  system  feedback 

configuration.,  The  problem  then  reduces  to  Case  1. 

In  the  above  manner  the  majority  of  feedback  problems  can  be  put 

into  zero  leakage  category.  For  those  situations  where  a  parallel 

path  from  the  overall  system  input  node  to  the  output  node  (from  I 

to  0)  exists,  the  following  results  are  obtained: 

Since  T  ■  t  .  +  kt  .t  1(1   -  kt  ) 

oi     ci  os/       cs 


using  Equation  (3=11)  and  notation  (3=12)  we  get 

,*    AT/Tf        i  /,   Jot 

\    •     Ak/k£      (1  -  kot  )   ^   "   t 


) 


fofcH 


sl  *  (  -ttV-  )  I1  -  -^-)  <3-18> 


Equation  (3-18)  can  also  be  written  as 

>t  J 

1  +  I 
o 

and  from  (3=19)  by  using  Equation  (3=17)  we  get 


■J   oi'  -     1  (3=19) 


or 


and 


T 
o 

oi 
oi 

(k  /k.  )  +  L 
of             o 

T 
f 

1  +  L 

o 

T 
f 

(1  +  L  )   -  (t    ./T  )(1  -  k  /k_) 

o                oi      O                  O      I 

T 
o 

L    +  k  /kf 
o         o 

At 

s 

T£  -  T                     (1   =  t    ./T  )    (1   - 
f         o                               ©i     o 

W 

T 
o 

T                                   L    +   (k   /k-) 
o                                 oof 

(3=20) 


III-3  -  Further  Investigation  into  the  Meaning  of  Sensitivity  Function, 

In  prev  ious  discussions  we  defined  the  system  sensitivity  as 

T         1 

(3=21) 
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where  L  is  defined  as  the  open  loop  transfer  function,  and  T  is  the 

closed  loop  transfer  function.  Let  us  now  investigate  the  value 

dT/T 
dL/L 

Referring  to  Figure  3-5  we  have 

G        1    GH 


1  +  GH       H  1  +  GH 
But  L  *  GH     and       dL  -  HdG 


dT/T    -  _L  dT_   -  _L  __dT_ 
dL/L       T  dL       HT   dH 


1  +  GH       1  +  L 

Now  by  definition  dT/T   D  _C  i.e.,  sensitivity  of  the  closed-loop 

dL/L      b0 
to  open-loop  perturbations,  and  comparing  Equations  (3-21)  and  (3-22) 

Weget  S*   -  si  (3-23) 


Also  from  relation 


S„  ■   Si 


'k    °0      i  +  l 

we  can  think  of  the  sensitivity  functions  as  the  error  transfer  func- 
tion for  a  unity  feedback  system  with  forward  elements  GH,  because 

R  -  C  ■  E 

or  RGH 

R  "  1  +  GH   "  E 

then  E        x  x 


R      1  +  GH       1  +  L 
In  control  system,  the  error  transfer  function  is  always  desired 
small,  ideally  zero,  and  we  may  therefore  observe  that  a  small  sensi- 
tivity function  is  desirable. 

Since  the  sensitivity  function  is  dependent  upon  frequency  and 
is  a  complex  number  for  any  arbitrary  frequency  we  must  further  inquire 
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into  its  meaning. 

Let 

then 


or  sc 


*  1 

(1  +  Pcos  0)  +  j(  Psin  0) 

_T  1  +  Pcos  0  ,  T    sin  | 

m      S  m  .    j  -E - 

i  +  aPcos  ^  +  r  i  +  aTcos  0  +r^ 

(3=24) 


or 


,c  eT  1 

>  »        J)  a  'i 

o  K  I " T"?~  °j 


ysn°         fsin  0 
1  +  2Pcos  0  +  T  "  1  +  P  cos  0  J 

(3-25) 
It   is  desirable  to  interpret  both  of  the  above  equations   from 

a  physical  view  point.     We  then  proceed  in  this  direction. 

c 

A  -  A  Synthetic  Interpretation  of  S 
•  o 


We  define  the  following  four  quantities: 

,M  D   d|T|/|T| 
*M  "   dJLJ/JL 


per  unit  increment  in  I  T|       /*.9M 
"LI       per  unit  increment  in  j  L  \ 

A  D   d  ^  /  T 

S.  e ■   increment  in  *— —  (3-27) 

d  J_L  increment  in  [   L 

M  D   dlTI/lTl       _£®r  unit  increment  in  I  T  j_      ,3  2g* 
A  "     d  /JL  increment  in  /  L  IJ-/OJ 

S   ■    d  *-—  *       increment  in  ^—  (3-29) 

dlL)/lL|       per  unit  increment  in  |  L  | 

Where  |'a|  and  _/_a  indicate  the  magnitude  and  phase  angle  of  the 

complex  number  a,  respectively.   It  is  necessary  to  compute  I  T  j 

and  (T     in  terms  of  |  L|  and  /L  in  order  to  evaluate  these  functions 

L  -  Pej^ 
Therefore  |  L  |  «  P     and  £L_    -  0  (3-30) 
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and 


H       1  +  ?ei9 


let  „    .     |H|eJ>!B- 

therefore 

T  ■ 


and  G     -     Ide^0- 


(3-31) 
(3-32) 


P 


lH,v/l  +  2Pco.  ^  +  r2 


^[fi-tan"1       Prtm#        -  Zg] 

e  u  i  +rcos  9        -J 


J 

(3-33) 


We  are  now  in  a  position  to  compute  the  functions  defined  by  Equations 

(3-26)  to  (3-29).  Note  that 

P      -    |  Hi   •  \G|  dP      -    |  H|  «d|G|      (3-34) 

0      */H_  +  2<L  dd    -    d  /G  (3-35) 

From  Equation  (3-26)  ve  have 


.M  dlTl/lTl 


M 


dlL|/|L\ 


-    I  nl    I I 15"       dL  lHl\/n»  2Tcos  d+P2""  , 


P 


(3-36) 


But  | H  |  Is  not  a  function  of  the  variable  part  of  P  then 

s*  -  \/i  +  2  r  cos  2  +r2 

H 

1  +    Pegs  0 
1  +  2TCOS  0  +  P2 

By  similar  technique  we  find 


i    path    vi     ■     iucu 

r      r       1 

■  n/i  +  2rCos  0  +  P2J 


(3-37) 


1+P 


C08 


'A 


,M 


1  +  2Pcos  0  +  P2 


A       i  +  2rcos0  +  r2 

-  I  sin  j 


<- 


1  +  2Pcos  ^+  P' 


(3-38) 
(3-39) 
(3-40) 
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1<1 


s      a*,     ■ 


<jM     _M     -A 

V    V   V 


ideally  sesr©.     Tc 

Che  mapping  of  the  loci  © 

(hence  all  four)  will  be 

B  -  Mapping  ©f  S. 


tniiusBiesfjLcs 


ar  plane 


We  require 


Where  M  is  a 


+  Pcos  4 


4=  2  P  cos  $  +  P 


1  + 


or 


t 


1  +  2se  +  3S  + 


+  (1 


Thus  the  loci  of  constant  §1  are 

coordinates  1 

x      2M 

and 

A 
C  -  Mapping  of  S  Loci  onto 


We  require 


where  Q  is  a  constant. 


|_  sin  0 


I  +  2  P  cos  0  +  P 


2         1 
+  x     +  y 


or 


+  1)      + 


/JL-\ 


(3-43) 


(3=44) 


(3-45) 


(3=46) 
(3-47) 


(3-49) 


(3-50) 
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Thus  the  loci  of  constant  S  are  a  family  of  circles  with  center 

coordinates  x  ■  •  1,   y  »  -  1        (3-51) 

2Q 

and  radius  -  _1_  (3-52) 

2Q 

Both  families  of  circles  derived  for  sj*  and  S  are  plotted  in  Figure 
3-7  and  it  is  observed  that  the  region  of  poorest  sensitivity  is 
near  the  (-1,  0)  point,  the  sensitivity  improves  in  every  radial 
direction  from  that  point. 
D  -  Mapping  of  J  S  ]  Loci  onto  the  Polar  Plane. 


From  Equation  (3-25)  we  require 
1 

7l+  2  Peo»  /i  +  r2 


M  (3-53) 


where  M  is  a  constant.  In  cartesian  coordinate  it  is 

(x  +  l)2  +  y2  -  /  1  \  2  (3-54) 


■(+) 


The  loci  of  constant  Sc  lare  a  family  of  circles  with  center 

coordinates,  x  »  -1,  y  ■  0  (3-55) 

and  radius  ■   1 n 

M 

Sc 
E  -  Mapping  of  /  o  Loci  onto  the  Polar  Plane. 

From  Equation  (3-25)  we  require  that 

tan-l     r«l«  L  .  c«  (3-57) 

1  +  p  cos  0 

or  in  cartesian  coordinate 

y  -  tan  C  (1  +  x)  -  C(l  +  x)        (3-58) 

The  loci  of  constant  phase  angle  of  the  sensitivity  function  are 

therefore  straight  lines  through  the  point  (-1,  0)  with  slope  C. 

This  family  of  lines  comprises  the  orthogonal  trajectories  to  the  loci 

of  constant  I  S 
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F  -  Result  of  the  Interpretation. 

For  any  point  on  the  polar  plane 


(<y  *  w 


so 


(3-59) 


where  each  term  is  evaluated  at  the  point  in  question.  Thus  the 

M      A 

magnitude  sensitivity  function  is  a  vector  measure  of  S  and  S  at 

every  point  of  the  polar  plane. 

The  only  extent  to  which  the  phase  angle  loci  will  he  of  use 

comes  from  the  consideration  of  the  phase  angles  zero  and  "]!_  rad. 

2C 
Along  each  of  these  lines  the  Imaginary  and  real  parts  of  SQ  respec- 
tively are  zero.  We  conclude  that  along  the  phase  angle  locus  corr- 
esponding to  zero  rad,,  magnitude  changes  in  L  cause  only  phase  changes 

in  T.  The  reverse  is  true  along  the  phase  angle ~n   rad.  locus. 

2 
Since  frequency  response  design  work  is  generally  carried  out  on 

the  amplitude-phase  plane,  the  mapping  of  the  magnitude  sensitivity 

function  loci  onto  this  space  will  now  he  considered. 

III-4  -  Mapping  of  the  Magnitude  Sensitivity  Function  Loci  onto  the 

Amplitude-Phase  Plane. 

The  magnitude  sensitivity  function  loci  can  be  mapped  onto  the 
amplitude-phase  plane  simply  by  inserting  the  variable  change 

Z  ~JU? 
into  Equation  (3-53) .   There  results  from  this  substitution, 

e2Z  +  2eZ  cos  0  +  1  -  t   1  \   2       (3-61) 
Plotting  this  on  the  Z~$   plane  gives  the  desired  mapping. 

The  ultimate  purpose  for  the  creation  of  this  mapping  is  to  use 
it  as  an  overlay  in  design,  just  as  is  used  in  the  Nichols  Chart.  By 
a  simple  manipulation  it  will  be  shown  that  the  Nichols  Chart  Itself 
is  already  calibrated  for  our  use.   In  previous  sections  we  concluded 
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<- 


■K  -0      1  +  L  <3-«> 

where,  L(s)  •  G(s)  H(s),  Figure  3-5.  It  is  intended  to  point  out 

that  this  sensitivity  function  is  obtainable  directly  from  a  convent- 
ional Nichols  Chart,  given  the  open-loop  frequency  locus 

L(jto)  -  G(j6>)  H(jtl) 
for  the  linear  system  in  question.  The  procedures  is  quite  direct. 
Both  numerator  and  denominator  of  Equation  (3-62)  may  be  divided  by 

L(s)  giving, 

(3-63) 


S 


ST  - 
SK 


JL£si 


-1 


1  +  Us) 

We  seek  the  polar  coordinate  mapping  of  lines  of  constant  magnitude 
sensitivity  function  given  by 


K 


-  M 


or 


±L*1 


-l 


1  +  L(s) 


-1 


M 


let 
therefore 


L(s) 


-1 


-  s  +  jy 


x+  jy 


1  +  x  +  jy 


and  from  Equation  (3-66)  ve  get 

.2       \  2 


x  + 


M 


M2-l 


+   y    . 


-   M 


M 


M2-l 


(3-64) 
(3-65) 

(3-66) 
(3-67) 


We  recall  the  equation  for  the  M  circles  of  conventional  design 
(ref.  9)  and  note  that  it  is  identical  to  Equation  (3-67).  Thus  we 
have  the  interesting  result  that  the  M-circles  are  to  the  open-loop 
frequency  response  curves  as  the  magnitude  sensitivity  function  loci 
are  to  the  inverse  open- loop  frequency  response.  The  Nichols  Chart 
is  just  a  plot  of  the  M-circles  (as  well  as  the  N-circles  of  constant 
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closed -loop  phase)  on  the  amplitude-phase  plane,  therefore  the  mapping 
of  the  magnitude  sensitivity  function  loci  onto  the  maplitude^phase 
plane  is  the  Nichols  Chart,  provided  that  we  read  intersections  as 
derived  from  inverse  open- loop  frequency  response  curves. 

The  combination  of  the  logarithmic  amplitude  scale  and  linear 
phase  scale  on  the  amplitude  phase  plane,  with  the  symmetric  Nichols 
chart  is  shown  on  Figure  3-9.  Consider  a  typical  open-loop  locus  and 

its  inverse,  on  the  extended  Nichols  chs^S  as  in  Figure  3=9,  The 

»1 
inverse  function  L  (jto)  is  obtained  @©  &   point  by  point  basis  graph- 
ically and  is  symmetric  with  the  direct  function  L(jfo)  with  respect 
to  the  point  (0  DB,  o  deg) ,     The  Nichols  chart  is  shown  over  a  §40 
degree  range  and  is  seen  to  be  symmetric  (Ref.  9)  about  the  axis 
phase-0  deg. 

One  further  simplification  may  be  made  as  shown  in  Figure  3=10 

=1 

If  we  imagine  that  the  phase  calibration  for  L  (jb)  locus  on  the 

Nichols  Chart  is  opposite  to  that  commonly  used  then  the  same  M° 
contrours  used  in  ordinary  control  system  design,  measure  sensitivity 
function  magnitude  for  use  to  the  extent  of  being  properly  labelled. 
In  this  case  L(jO)  and  L~  (jw)  are  symmetric  with  respect  to  zero  DB 
line,  with  phase  measured  in  opposite  directions  for  each.  For  the 
particular  example  of  Figure  3=10  we  see  that  the  maximum  magnitude 
sensitivity  function  value  is  M  ■  2.0, 

If  it  was  desired,  the  magnitude  sensitivity  function  for  the 
closed  loop  system  as  a  function  of  frequency  could  be  read  directly 
off  the  Nichols  Chart  as  the  intersections  of  the  M-contours  with 
the  inverse  open- loop  frequency  response  curve. 
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The  simplification  introduced  by  foregoing  interpretation  of 
the  sensitivity  function  magnitude  and  phase  loci  could  be  very  use- 
ful In  the  design  of  linear  systems  with  regards  to  sensitivity 
specifications. 
Special  Case: 

If  the  forward  elements  can  be  separated  into  two  tandem  parts, 
one  of  which  displays  sensitivity  to  the  environment  (denoted  by  an 
"environmental  parameter"  x),  viz., 

G(s,  x)  -  G^s)  G2(s,  x)  (3-68) 

Then  we  may  compute  the  sensitivity  of  the  closed  loop  locus  of  the 

T 
overall  system  to  the  environmental  parameter  (denoted  S  )  according 

to 

C      dT/T        dT/T   dx/x 
0  "   dL/L        dx/x   dL/L 

-  S*  S*  (3-69) 

£  .  £  -L-  .  Sjj  S°  (3-70) 

o 
In  the  present  example 

dL/L      dG2(s,x)/G2(s,x) 

5x  "   dx/x   "        dx/x 

-  S  Z  (3-71) 

x 

Therefore        _,     r   G, 

ST  -  S^  S  l  (3-72) 

x      0   x 

where  the  first  term  of  the  product  is  evaluated  graphically,  and 

the  second  term  is  found  by  differentiation.   It  should  be  noted  that 

if  the  parameter  x  appears  in  G2(s)  as  GAa)     -  xG  (s),  then 

G2  T     c 

S  '      ■  1  and  S   ■  S_  ,  which  is  in  accordance  with  our  relation 
x  x     0  ' 

C      T 
S   -  S    where  k  appears  as  in  Figure  Al-1  of  Appendix  1. 
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Fig.  3-10:   L  and  L*  on  the  Nichols  Chart. 
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III-5  -  Philosophy  of  the  Frequency  Response  Approach  to  Sensitivity 
Problem. 


In  previous  sections  we  discussed  how  sensitivity  and  frequency 
are  related  and  how  by  utilising  the  Nichols  Chart,  the  magnitude  and 
phase  of  the  sensitivity  function  and  its  peak  magnitude  can  be  obtained 
Another  approach  which  could  be  of  high  importance  in  most  design  cases 
is  discussed  as  follows: 

Sensitivity  specifications  could  be  given  in  terms  of  the  systems 

behavior  on  the  jto  axis.   It  is  possible  that  the  allowed  variation 

in  T(jlO)  may  be  exceedingly  small  over  one  frequency  range  and  fairly 

large  in  other  frequency  ranges.   The  specification  may  obviously  be 

in  a  wide  variety  of  forms.   For  example  the  nominal  T  ,  may  be  as 

shown  in  Figure  3-11,  with  bandwith  U). ,  amplitude  peaking  M  ,  and  the 

l  o 

specification  may  dictate  that  despite  parameter  variations,  the 
magnitude  peaking  should  never  exceed  30  percent,  the  bandwidth  should 
be  in  the  range  0,8(o,  to  1.2W  ,  and  that  from  zero  frequency  to 
0.5CX,  )t|  should  never  be  less  than  0.95.   Finally,  some  thought 
should  be  given  to  the  intermediate  and  the  far  off  frequency  ranges. 
For  example,  the  behavior  shown  by  dotted  lines  in  Figure  3=11  is 
usually  intolerable,  as  it  results  in  high  frequency  ringing.   A 
reasonable  statement  might  be  that  for  1.2&  <  6><  ft.  ,  |T(j6))|<  3DB 
and  for  (J  >  U).,   T(jW)   must  not  exceed  a  predetermined  value. 
The  basic  philosophy  of  the  frequency  response  approach  to 
sensitivity  reduction  is  deduced  from  Equation  3°17,  rewritten  here 
T         (k  /k.)  +  L 

O  O   I o  /<1  =,«» 

-        m  (3~73) 

X  o 

This  equation  applies  whenever  the  leakage  transmission  t   is  zero 
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0 

0  DB. 

1 

W3 

/ 
/ 

-3DB. 

V 

-V  -     i- 

V      --  i 

1           \     V-.      -o" 

Fig.  3-11:   Undesirable  High 
Frequency  Peaking  in  T(jk>)  to 
Parameter  Variation. 


-L  (jU»)  -  A  +  JB 
o   x 


Range  of  P  (ju>) 
oVJ  x 


Fig.  3-12:   Showing 
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Fig.  3-13:   To  Obtain  Permissible  Location  of  -L  (\Gi  ) 

o  J  xy 


Range  of  P  /P(ju>  ). 
o     x 


77 


and  therefore  could  be  used  in  most  control  problems.   Suppose  that 

at  some  specific  frequency  b  ■  j(J,  plant  parameter  variations  are 

such  that  (k  k.)  at  8  °  j(J  may  lie  anywhere  inside  the  region 

indicated  in  Figure  3-12.   Suppose  also  that  -  L  (j^  )  is  given  by 

the  complex  number  A  +  jB  located  at  Q.   Then  Equation  (3-73)  states 

that  _ 

JjL.    m     _9JL 
f       QN 

The  range  of  variation  of  the  vector  QV  fixes  the  range  of  variation 

of  T  /T„  at  s  »  jo  .   Suppose  for  example,  that  it  is  required  that: 

1.20  >  T  /Tr  >  0.80 
o  t 

at  s  ■  jt)  .  '  The  problem  is  then  to  find  the  locus  of  -L  <j#  )  which 


barely  satisfies  this  inequality.   A  suitable  procedure  is  as  follows: 

>  0.80  with 


Suppose  we  seek  the  locus  of  the  boundary  of I   o 


Tf 
point  N  ■  1  +  j  as  center,  draw  a  circle  G,  of  any  radius  R  (Figure 

3-13).  Next  draw  a  circle  C.  of  radius  0.8R,  using  as  center  the 
point  (the  boundary  of  k  /k.)  which  appears  to  be  the  closest  to  the 
first  circle  -  this  appears  to  be  point  J  in  Figure  3-13.  The  inter- 
section of  the  two  circles  determines  a  point  Q.   As  a  final  check, 
use  Q  as  center,  and  check  whether  a  circle  of  radius  0.8R  cuts  k  /k^ 
locus  only  at  J.   If  so  Q  is  on  the  boundary  of  the  locus  of 
T  /Tf  /  >  0.80  at  s  ■»  j£>  .   The  above  construction  is  repeated,  using 
different  values  for  radius  R,  until  the  locus  of  Q  is  obtained.   Such 
a  locus  for  1.20  >  T  /Tf  I  >  0.80  is  shown  in  Figure  3-13.   Clearly, 
"L  (jt)  )  must  be  located  inside  the  indicated  region  in  Figure  3-13. 
On  the  other  hand  if  there  is  the  additional  requirement  that  the 
phase  of  T  /T,  must  not  change  more  than  20  ,  then  the  permissible 
range  of  -L  (jO)  is  reduced  to  the  smaller  region  shown  in  Figure  3-13, 
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It  is  extremely  important  to  note  that  due  to  Equation  (3-17)  the 
sensitivity  requirements  essentially  determine  the  gain  and  band- 
width  of  L  (s)  with  hardly  any  need  to  refer  to  the  specific  feed- 
back  configuration  that  is  used.  At  best  one  configuration  may  be 
more  efficient  than  another  in  the  practical  realisation  of  the 
required  L  (s).  The  latter  is  actually  a  problem  in  network  synthes: 


In  case  1  we  assumed  that  in  the  fundamental  feedback  equation 

T  »     t.+kt.t  /(l  =  kt  ).  the  leakage  transmission  t  .  was 
oi     ci  os       cs  °  oi 

negligible  in  comparison  with  T,  over  the  range  of  variation  of  k. 
As  previously  discussed  this  is  not  always  the  case,  When  the  leak- 
age transmission  t  „  was  not  small  in  comparison  with  T,  then  Equations 
(3=18)  and  (3-20)  were  resulted  which  are  rewritten  here; 


Tf 


L  )  (3-74) 


(3-75) 


(3-76) 


T 
o 

where  Tf  is  the  value  of  T  when  k  has  the  value  k  »  and 

AT     *     T-  -  T 

f         o 

Ak     s     k„  -  k 

f         © 

Now  t  .   is  a  transfer  function  which  is  independent  of  k.   If  t  '  is 
known  and  the  desired  T  is  known,  then  one  may  relate  the  maximum 
permissible  range  of  variation  of  T9  over  any  frequency  range,  to  a 
corresponding  maximum  range  of  variation  of  T  -  t  - .  In  other  words 
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(a) 


tol/T  freq.  range 


-tQl/To  (1  -  ko/kf)  freq.  locus 


o 


locus  k  /k- 
o  f 


1  -  'ol/To  (1  -  ko/k£> 


Fig.  3-14:   Construction  for  finding  T  /T_,  when  t  .  is  not  zero. 

o  t       oi 
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one  can  work  graphically  with  Equation  (3°75)  in  exactly  the  same 

manner  as  Equation  (3-73)  was  used  in  Case  1..  The  only  difference 

is  that  one  must  first  decide  what  is  the  permitted  range  of  vari° 

ations  of  T,  =  t  .,  rather  than  that  of  T„. 
f    oi  f 

The  above  is  satisfactory  if  t  ,  is  known ,  On  the  other  hand, 

©i 

there  may  be  situations  when  t   is  not  precisely  known  until  part 
of  the  design  has  been  done.   In  such  a  case  some  cut  and  try  is 
inevitable.  The  following  modification  of  the  graphical  procedure, 
suggested  by  Horowitz  is  an  alternative  method  of  handling  the  non- 
zero leakage  transmission  problem.  Equations  (3=75)  to  (3=77)  are 
used.   Sketch  the  frequency  locus  of  t  /!  (Figure  3°14a)  in  the  com« 
plex  plane,  over  the  frequency  range  in  which  control  of  variation 

in  T  is  desired.  Also  sketch  the  loci  of  k  /k,.  From  these  two 

©  f 

figures,  obtain  the  loci  of  =(1  -  k  fkj   (t  ./T  )  as  in  Figure  3-14b. 

O   £     Oi   O 

Finally  combine  the  latter  and  the  sketch  of  k  /k£   into  one  figure 

o  t 

(Figure  3°14) .  In  Figure  3°14e  we  have 

T 

_j_         m 

T    S   MB 
o 

in  accordance  with  Equation  (3=76),  The  design  problem  is  to  select 

the  loop  transmission  L  so  as  to  satisfy  the  specification  of  T  /T_. 

III-6  -  Use  of  Polar  Plot  in 

■! 


The  relation,      „T 


0      1  +  L  (s) 

© 

suggests  that  the  inverse  polar  plot  could  be  useful  in  studying  the 

sensitivity  function.  Above  relation  could  be  written  as 

=1 

i  -  *$  ■  ~rhr  m     L°  -I  (3~78> 

K  °  1  +  \  1  +  L       (s) 
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Referring  to  Figure  3-15  it  is  noticed  that  if  OF  represents  L  "  (j&>), 

then  BP  -  1  +  L  "  (jG>)   and,  x 

C  Lo"   (jW)  OF  j<* 

so   -   — 2~n BPMc  (3"79> 

0      1  +  L   (j0>)    Br 
As  discussed  in  Equations  (3-63)  to  (3-67)  we  have 


"[t^tj  2  +y2 "  (ttt) 


(3-80) 


where       x  ■  real  part  of  L   (jto), 

-1 
y  -  imaginary  part  of  L   ( jtt) . 

M 

Equation  (3-80)  represents  a  circle  having  radius  ■  5—  ,  and 

1  -  M 

2 
center  at   -  M     on  the  x  axis.  (3-81) 

M2-l 

The  loci  of  constant  M  described  above  is  useful  for  determining  the 

value  of  M    for  a  control  system  for  which  L  (s)  and  eventually 

L   (s)  is  known.  The  radius  of  the  M  circle  and  the  distance  from 
o  v  ' 

the  origin  to  the  center  of  a  given  M  circle  are  both  a  function  of 
M  alone.  Referring  to  Figure  3-16,  where  P  is  the  point  of  tangency 
from  the  origin  to  the  M  circle,  one  observes  that  the  ratio  of  BP 
to  OB  is  a  constant  for  a  given  M. 

BP  -  radius  of  M  circle  ■  M/(M2  -  1), 


so 


OB  -  M2/(l  -  M2), 

M 

■ 

1 

M 

ci_  w     _  . 

1  -  M2 

8  in  t  ■ 

M2 

1  -  M2 

'  is  solely  defined  in  terms  of  M. 

OC       OP  cost 
OB        OB 

2 

OB  cos  Hf 

OB 

■ 

1  - 

1 
M2" 

M2-l 
M2 
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Complex  Plane 


+  Lo  «  >  fyVo  > 


Fig.  3-15:   Conventional  Form  of  Complex  Plane  Diagram  for 
Control  System  of  Equation  (3-78). 


Fig.  3-16:   Geometry  of  M  Circlet, 
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2   2 
or       OB/OC  -  M  /(M  -  1). 

Hence  C  represents  the  location  of  -1  +  jO,  point  to  the  same  scale 

2 
for  which  BP  represents  the  length  M/M  -  1. 

For  a  given  L   (s)  function  to  have  a  maximum  M  value  of  M  it 
o  m 

must  also  be  tangent  to  the  M  circle  having  a  value  of  M  .  By  using 

m 

these  facts  the  following  procedure  can  be  formulated. 

1  -  Draw  the  locus  of  L   (jo)  function  on  the  complex  plane 

with  the  value  of  K  -  1,  K  being  the  gain. 

2  -  With  the  value  of  M  specified,  determine  the  angle  *f  from 

m 

S.  -   -  "l 


-(+) 


3  -  Using  the  angle  ^  determined  from  2,  draw  the  line  OD  on  the 

complex  plane. 

4  -  By  cut-and-try,  draw  the  circle  that  has  its  center  on  the 

negative  real  axis  and  is  tangent  to  both  L   (jo)  locus 

o 

and  to  the  line  OD  at  some  point  P  as  shown. 

5  -  Draw  a  perpendicular  PC  from  the  point  of  tangency  P  to 

the  negative  real  axis. 

0C(1/K)  m   1.0 

then  K  -  OC 

This  value  of  K  is  the  gain  associated  with  L  (s)  that  will  yield  the 

specification  of  M  .  Knowing  the  proper  net  gain  K,  one  may  then 
m 

change  by  K  the  scale  of  the  plot  used  in  preceding  steps  and  evaluate 
the  sensitivity. 
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IV 
SPECIAL  TOPICS 
IV«=1  -  Discussion, 

In  previous  sections  we  discussed  the  sensitivity  function  and 
shoved  that: 

w—  •* 

T         1  *- 

(4-1) 


t[ 


<-TT    ■--* 


where  F,  is  the  return  difference  with  respect  to  the  parameter  k. 

The  parameter  k  could  be  any  parameter  such  as  gain  K,  pole  P  ,  or 

zero  Z  ,  and  t   is  the  leakage  or  direct  transmission  between  input 
1       ox 

and  output ,  as  shown  in  Figure  4=2.  Although  the  sensitivity  function 
has  been  defined  and  several  of  its  properties  examined,  it  is  not  yet 
apparent  how  one  may  utilise  it  in  system  design,  or  more  fundament- 
ally on  what  basis  its  characteristics  should  be  chosen.   Similarly 
in  order  to  compare  the  relative  sensitivles  of  different  parameters 
it  is  convenient  to  have  an  overall  measure  of  sensitivity  which  is 
independent  of  either  time  or  frequency. 

IV-2  -  Specification  of  the  Sensitivity  Function. 

A  common  design  approach  consists  of  making  the  "loop  gain"  or 
return  difference,  extremely  large  in  the  band  of  frequencies  in 
which  the  input  signals  are  expected  to  be.  Then  the  open  loop 
response  is  caused  by  appropriate  shaping  to  fall  off  as  rapidly 
as  possible  to  unity  gain  outside  this  band.  A  maximum  roll  off 
rate  of  about  33  db  per  decade  permits  unconditional  stability  to  be 
achieved,  while  a  greater  rate  results  in  conditional  stability.  Thus 
realization  of  low  response  sensitivity  by  making  the  magnitude  of 
the  sensitivity  function  small  in  the  desired  band  of  frequencies 
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input 
R(s) 


Fig.  4-1:   A  Basic  Feedback  System. 


Output 
C(s) 


(a) 


-(b(8)[M(s)/K]) 


R(s)# 


A(s)B(s)[M(s)/K] 


C(s) 


(b) 

A(s)B(s)fM(s)/(S  +  Zj)l 
(1  -  B(s)[M(s)/(S  -  Z 


-B(s)fM(s)/(S  ±   Zl» 
1  +  B(s)[M(s)/(S  +  Zt)] 


Us) 

SA(8)B(s)[M(s)/(S  +  Zt) 
1  +  SB(8)[M(s)/(S  +  P1)J 

Fig.  4-2:  Alternate  System  Configuration  of  Figure  4-1. 
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(c) 


Af«)R(fi)[M(fi)  <.S  -I- 


[S  +  B<8)]  [M(«)  (S  +  P^] 


R(s) 


^ 

S  +  B(s)[M(s)  (S  +  P^] 


A(s)B(s)[M(s)  (S  +  P^] 
S  +  B(a)  (M(s)  (S  +  P^] 


C(s) 


(d) 


Input 
R(s) 


Out? 
C(s 


Fig.  4-2:   Continued, 
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inevitably  results  in  the  problem  of  system  stability.   In  the 

following  several  criterions  for  the  specification  of  sensitivity 

functions  are  discussed. 

(a)  -  The  Forced  Response  Criterion. 

Suppose  the  desired  response  in  the  time  domain  is  given  by 

c(t)  and  the  response  error  by  e(t).   Referring  to  Figure  4-1,  if 

(s  +  Z  )  (s+  Z  )  -— 
M(s)  -  K  ,.  ;  p\  ,.  T  A  ..__  (4-2) 


(s  + 

pl> 

(8   + 

P2>    "" 

C(s) 

-  . 

A(8) 

B(s)  M(s) 

R(s) 

1  + 

B(s)  M(s) 

then  T<»>  -  -ifit  -  Ti  IK  mi  (4-3) 

and  assuming  variable  element  x  is  associated  with  M(s)  only,  then 

sT/.s  .    x      dT(s)   .  x     dM(s)/dx 

V8'      T(s)      dx       M(s)   1  +  B(s)  M(s) 

(4-4) 

for  x  »  K,  Z.  or  P.  Equation  (4-4)  becomes, 

(4-5) 

(4-6) 


SK(S>     " 

sz(8)    * 

J 

Sp(8)       « 
^i 

1 

1  +  B(s)  M(s) 

1 

zi 

1  +  B(s)  M(s) 

i 

s  +  P 

1 

+ 

B(»)  M(s) 

(4-7) 

The  basic  problem  is  determination  of  criteria  for  the  specifi- 
cation of  the  sensitivity  function  for  an  arbitrary  input.   After 
the  form  of  the  desired  transmission  function  has  been  ascertained, 
based  on  some  dynamic  response  criterion,  Equations  (4-5)  to  (4-7) 
reveal  that  the  system  parameters  which  remain  to  be  specified  are 
the  zeros  of  the  sensitivity  function.   The  poles  of  the  sensitivity 
function  and  the  transmission  function  are  identical,  except  for  the 
added  poles  Z  and  P  ,  provided  the  product  A(s)  B(s)  M(s)  in  Equation 
(4-3)  is  chosen  so  that  no  poles  are  added  to  those  of  the  reciprocal 
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of  F .  (s).   The  forced  response  criteria  determines  the  sero  of  S(s) 
so  that  e(t)  behaves  in  the  desired  manner  over  the  interval  or 
intervals  of  interest  in  the  time  domain.   A  considerable  literature 
exists  whereby  one  may  relate  time  domain  responses  to  the  complex 
frecuency  domain  characteristics  of  a  given  network  configuration. 
These  results  are  useful  when  the  inputs  are  restricted  types  and 
when  the  system  poles  are  simple.   Functions  of  third  or  higher  order 
with  multiple  poles  and  possible  poles  on  the  jw  axis  are  difficult 
to  deal  with   If  the  sensitivity  function  is  to  be  specified  only  on 
the  basis  of  tKe  forced  response  to  a  given  input,  then  the  zeros  of 
S(f )  are  chosen  as  follows- 
if  T  (s)  consists  of  a  number  of  inputs  of  the  form, 

then  the  system  output  is  given  by, 

Cj(8)  "  aj  sA   T(8)  (A°9) 

and  the  forced  response  is 

c  (t)  ■*  a^-s^e'V  (4=10) 

Taking  the  derivative  of  c.(t)  with  respect  to  a  variable  parameter 

-s.t   dT(-s.) 

dx  J  (A°U) 


x  we 

i  get, 

dc^t) 

■ 

dx 

V 

But 

from  Equation  (4-4) 

we 

have 

dT(-8j) 

■     s 

-  ST  (~s,)/x  (4-12) 

dx  r.  Y 

Thea  substituting  Eauation  (4-12)  into  (4-11)  we  get 

te  At)  a  -st 

1 .  .  _1-  t(-s  J  sVO  e  J        (4-13) 

dx         x        j   xx   j 
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T 
Thus,  if  the  sensitivity  function  S  (s)  has  a  zero  at  the  pole  of 

the  transform  of  the  input  function,  the  system  forced  response  is 
invariant  with  changes  in  x.   If  s  '  s  are  expected  to  lie  in  some 
known  range  of  values,  the  sensitivity  zeros  may  be  distributed 
accordingly. 

The  above  result  is  also  valid  when  s  is  imaginary  but  in  that 
case  both  the  steady  state  amplitude  and  phase  responses  are  of 
interest.  For  a  sinosoidal  input,  the  forced  response  and  its  deriv- 
atives with  respect  to  x  become 


c  (t) 


r  '      8 

r 


T(jsr)|  sin  srt  +  0  (jsr)    (4-14) 


where  .1  ImT(js  ) 

0(js  )  -  tan 


and 


r'  ReT(jsr) 


dcr(t}  *-        "  JT(jsJ| 


dx         8 


^  ^sin  (srt  +  f)  «■  dx 


d0(js J 


T( jsr)  £-^-  cos  (sj:  +  0)\  (4-15) 


Writing        T(js  )  - 


T(jsr) 


»rc  +  m 


(4-16) 


Then  T      .  x      df  It|  ej*] 

xu  r>        -TT3F         dx 


^TT   dx    +  JX  dx 

-  SR(jsr)  +  jSI(jsr)  (4-17) 

The  R  and  I  subscripts  in  the  last  equation  refer  to  the  real  and 

imaginary  parts  of  S(js  ),  respectively.   Thus  comparing  Equations 

(4-15)  and  (4-17),  if  zeros  of  S(s)  are  assigned  at  +  js  ,  then  both 

S„  and  S„  are  zero  and     r^  '     becomes  zero  and  is  invariant  with 
R      I  £— 

changes  in  x. 
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In  this  -action  we  demonstrated  how  variations  in  the  forced 

response  may  be  controlled  by  suitably  specifying  the  sensitivity 
function,   A  design  based  only  upon  the  forced  response  is  likely  to 
result  in  an  unsuitable  selection  of  the  zeros  of  the  sensitivity 
function,  if  the  system  is  of  high  order  and  the  input  is  of  a  known  in- 
variant form.   For  example  if  a  fifth  order  system  is  to  be  subjected 
to  a  sinosoidal  input,  then  five  sensitivity  zeros  must  be  assigned. 
Two  zeros  are  given  values  corresponding  to  the  poles  of  the  input 
transform,  there  is  no  reasonable  basis  for  selecting  the  remaining 
three.  Furthermore,  if  interest  is  centered  upon  the  part  of  the 
system  response  where  the  transient  has  an  appreciable  value,  then 
even  the  forced  response  criteria  outlined  have  little  value  in  deter- 
mining suitable  sensitivity  zeros.   Then  another  technique  based  on 
"minimum  mean  square  error"  was  proposed. 
(b)  -  The  Minimum  Mean  Square  Error  Criterion 

The  minimization  of  the  mean  square  difference  between  a  desired 
and  an  actual  response  is  fully  discussed  in  literature.  This  tech- 
nique is  directly  applied  by  Mazer  to  the  sensitivity  problems,  its 
main  advantage  is  that  the  problem  is  solved  entirely  in  the  complex 
frequency  domain,  thus  the  necessity  for  obtaining  the  error  in  the 
time  responses  of  hi'b  order  feedback  systems  in  terms  of  the  unknown 
zeros  of  the  sensitivity  function  is  avoided.   Before  proceeding  with 
details,  the  following  definitions  are  required: 

T(s)  -   gffi  (4=18) 


AT(s)  -  "-  T(s)  3(s)  (4-20) 
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E(s)   -  R(s)AT(s)  (4-21) 

Where  F(s)  is  the  laplace  transform  of  the  input  and  A  notation  is 
used  to  indicate  a  small  finite,  rather  than  differental,  change. 

The  mean  square  error  is  given  by 

/T 
e2  -  lim   -~-  )      e2(t)dt  (4-22) 

T-»oo   *l        -T 

Here  e(t)  is  the  inverse  transform  of  E(s).   The  square  of  the  error 
is  integrated  over  a  time  interval  of  length  2T,  its  average  is 
taken  and  the  time  interval  is  then  allowed  to  become  infinite.   It 
should  be  noted  here  that  T(s)  is  the  transfer  function  of  the  system, 
while  T  Is  the  time  interval.   As  shown  in  reference  (2), 

e2  =  ~~       J  AT(s)  AT(-s)  f   (s)ds        (A-23) 
lj    -jco  r 

Where  <J)  (s)  is  the  power  density  spectrum  of  the  input  signal.   Sub- 
stituting Equations  (4=19)  and  (4-20)  into  (3-23)  we  get 


i     fJ00/*   \2   Nn(s)NA("8) 


2 

To  minimize  the  mean  square  error  e  ,  one  takes  its  partial  deriva- 
tives with  respect  to  each  unknown  coefficient  of  Nft(s)  and  sets  the 
results  equal  to  zero.   This  yields  a  number  of  linear  algebraic 
equations  equal  to  the  number  of  unknown  coefficients.   Furthermore 
if  N  (s)  is  of  even  or  odd  degree,  ail  coefficients  of  odd  or  even 
powers  of  s  must  equal  zero.   In  order  to  obtain  a  suitable  density 
spectrum  for  use  in  Equation  (4-24),  aperiodic  inputs  such  as  steps 
and  ramps,  must  be  modified  by  being  repeated  periodically.   Although 
the  system  is  then  designed  on  the  basis  of  this  fictitious  input, 
the  periodicity  of  the  repeated  input  signal  determines  approximately 
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the  Interval  of  time  over  which  the  mean  square  error,  for  the  case 

of  an  aperiodic  input  is  minimized. 

Example 

Design  a  feedback  system  in  the  configuration  of  Figure  4-1,  with 
the  following  closed  loop  transfer  function 

5 


T(s) 


(4-25) 


(a  +  s  +  1)  (s  +  5) 
Since  the  system  is  of  the  third  order,  NA(s)  may  immediately  be 
written 


3  2 

N_ (a)     ■     s    +  a„s    +  a,s  +  a 
0  2  1  o 


(4-26) 
If  the  variable  element  is  the  forward  gain     ,    the  mean  square  error 

bee  >  es • 


e      ■ 


1  \   m( A  k\     T(a)T(- 

2flj        Z\a*\    K    9     D(   )D(' 


a) 


8) 


2  N2       2,  2 

(as    +  a  )   -  s   (s 


I    0r(s)dsl  (4- 


♦ .,)«] 


27) 


Since  e     is  the  sum  of  two  integrals,   one  may  solve   for  a     and  a 
separately         Taking  partial  derivatives  with  respect  to  these 
quantities, 


a.2 


2tfj 


<joB   r~K/        D(s)D(-s) 

2(a2s     +  aQ)   s  0r(s)ds 


and  similarly 
fla 


[■ 

1  P^/fllA2     T(s)T(-s) 

2Tlj       /        V    K/        D(a)D(°a) 

-JCO     x      / 

2(a2s2  +  aQ)  0r(s)ds 


(4-28) 

■     0 


>  a2  +  a  )  0r(s) 

Solving  Equation  (4-28)   and   (4-29)    for  a.  and  a     homogeneous   set  of 
equations  results,    the  solution  which  is 

a_     »     a       ®     0 


To  solve  for  a, , 


*e2       . 
*al 


L_     f5JO°(lAJL\2       li*l  Uz*)      r.282(S2+a  )lrf   (s)ds  »  0 


(4-30) 


and  from  that 


Jj"  ;<:>;<:;>    cs2)(-82)Ms)ds 


21U      ^4 CO  D(8)D(-8)  x       '*    "    "r 

ai      ■      " "Vic*  <4"31> 

1  i       r  J06 


/J~-fSM-   <.x-.)*r(.)d. 


2TTJ      ^        D<s)D< 

Thus  a  is  independent  of  AKS  provided  the  latter  is  small  enough  for 
the  first  order  approximation  ^T(s)  given  by  Equation  (4-20).  Further- 
more ^  (s)  being  positive,  a  is  always  positive  and  the  integrands 
are  written  as  conjugate  products  in  Equation  (4-31)  to  emphasize  this 
fact.  The  character  o£  the  input  remains  to  be  specified  so  that 
4   (s)  may  be  evaluated.  Suppose  the  input  is  a  sinusoid  of  frequency 
B   »  1/2  rad/sec.  and  that  the  mean  square  error  of  the  system  response, 
over  one  cycle  of  the  driving  function  is  to  be  minimized  for  small 
changes  in  k.  To  generate  an  appropriate  repetitive  input  signal  let 
the  input  sinusoid  be  multiplied  by  a  square  wave  of  frequency 
0i  ■  1/4  rad/sec  and  unity  height.  The  Fourier  series  and  power 
spectrum  of  the  square  wave  are 

r  (t)  -  ~sr   <sin  ft  +  i  sin  |t  +  \   sin  yt  +  —  ) 
sqv  '  7T43     45     4 


and 


*sq(<0)  "  \W   [[«*-*>+*<  <*+*>] 


+  \\l^~  |)  +  5<»+  f)l 


+  ^[$((0- f)  +  S«o+  -§->] 

+  .........  (4-32) 
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Here  5(63  +  nW)  is  the  Dirac  delta  function.  If  this  square  wave 
—   o 

is  driven  through  a  fictitious  transfer  function 

F(s)  -    2  S  2   -    ,  S  ;  t  o-  (4-33 


s  +  / 


+(D' 


then  the  power  spectrum  of  the  input  is  given  by 

2 

,2 


[to  -  /2>  '     q 


With  0  (w)  given  in  Equation  (4=32),  the  intergrals  in  Equation  (4-31) 
sq 

are  evaluated  by  summing  twice  the  values  of  the  integrands  at  W  ■  —  , 

3   5 

t  ,  t  i  etc.,  since  they  are  even  functions  and  thus  always  positive. 

The  summations  are  continued  until  the  desired  accuracy  is  achieved. 

The  values  of  the  integrands  fall  rapidly  with  increasing  Cd  ,   since 

the  degrees  of  the  denominators  are  much  larger  than  those  of  the 

numerators.  Using  only  the  first  five  terms  in  the  Fourier  series  of 

Equation  (4-32)  and  substituting  T(s),  D(s)  and  jf   (s)  we  get 


1 
2HJ 


i^^C^T     <s>Va)\(s)ds  -  0.751(10^2) 

"J  o» 


(4-35) 


Irtir^   (s)2(-s)2^  (s)ds  -  i.238(io  ) 

D(s)D(-s)    *  '  *  '   rr     '  v    ' 

~Jas  (4-36) 


(s)(~s)$  (s)ds  »  2.511(10°2) 


2TTj    <     D(s)D(-s)   xw/x  'Tx 
-Jco 

Then  by  Equation  (4-31)  we  get 

a       -     0,404 

and  the  mean  square  error  is 


(4-37) 
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±_       yJ7AS\2       T(i)t(-)        .,,2^2       .2  a  (   )d 
ST1J        ^   (^  K^  D(a)D(=»)        lC   S   ><8+V    jPr<8>d8 


(4-38) 
Equation  (4-38)   cam  be  written  as 


2 

'       *     "2TTJ 


1  rj°°  Ak  V       T(s)T(~s)      /     3.  .  3.  j  ,   N. 


+       2TTJ 


2TTJ 

and  utilizing  relations  (4-35)  to  (4°37)  we  get 

e2  -  0.00149  l-^Yj  (4-40) 

The  sensitvity  function  is  then 


D(s)      («Z  +  s  +  1)  (s  +  5) 

The  above  techniques  show  how  complicated  the  computation  could 
be  and  in  many  cases  the  use  of  computer  is  the  quickest  way.  The 
purpose  of  the  above  discussion  is  not  a  complete  design  technique, 
but  to  show  how  different  approaches  are  made  to  the  problem. 

for  more  insight  into  this  problem  reference  2  and  11  are 
recommended „ 

IV-3  -  Equivalent  Isaput  P®irturliation  Signal 

As  discussed  in  Appendix  1,  for  any  linear  system,  the  transfer 
function  T(s,x),  relating  a  response  to  an  input  function  can  be 
expressed  in  the  following  bilinear  form  in  terms  of  a  real  parameter 
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x  of  the  system  and  the  polynomials  A(s)4  B(s),  C(s),  and  D(s) 


T(s,x) 


P(s,x) 
Q(s,x) 


»  +  xB(s) 


C(s)  +  xD(s) 

where  P(s,x)  ■  a(s)  +  xB(s) 

Q(s,x)  s  C{s)  +  xB(s) 
From  Equations  (4-43)and  (4-44)  we  have 

A(a)  ■  P{s,o) 
B(s)  -  —  (P(s,x)  -  P(s,o)] 
C(s)  ■  Q(8,o) 


D(s)  -  —  [Q(s,x)  -  Q(s,o)] 


(4-42) 

(4-43) 
(4-44) 

(4-45) 

(4-46) 
(4-47) 
(4-48) 


With  the  aid  of  these  equations  the  logarithmic  sensitivity  function 
may  be  expressed  directly  in  terms  of  the  polynomials  of  the  original 
transfer  function. 


S(s,x)  ■  S 


T.  frg»T(s,x) 


But 


and 


T(s,x) 


x       ^T(sBx] 
T(s,x)      dx 

x  Q(sflx) 


>»*> 


"b  »         bx 


A(s)  +  xB(s)  |(4 
[_C(s)  +  xD(s) 


(4-49) 
(4-50) 
51) 


[8)C(s)  -  D(s)A(s) 
[Q(s,x)]2 

Substituting  Equations  (4-50)  and  (4-51)  into  (4-49)  one  gets 


S(8,x)  ■  8   -  x 
x  »  '      x 


B(s)C(s)  -  A(s)D(s) 
P(sBx)  Q(ssx) 


] 


|s)Q(s,o) 
P(s,x)Q(s,x) 


P(s8x)Q(88x) 


TP(s.x)  -  P(s,o)1  Q(s,o)  -  FQ(s,x)1 
P(s,x)Q(s0x) 


[Q(s,o)       P(sto)  I 
Q(sBx)       P(s,x) 


(4-52) 
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Now  if  the  variable  x  changes  by  an  incremental  amount  x,  then  an 
approximation  given  by  the  first  terms  of  a  Taylor  series,  T(s,x), 
changes  by  an  incremental  amount  given  by 

AT(s,x)  ■  T(ssx  +  Ax)  -  T(s8x) 
!snx} 


-  (-*§-)  T(s,x)S(«,x)  (4-53) 

For  a  given  input  function  R(s)9  the  system  normal  response  is 
given  by  Y(s9x)  ■  R(s)  T(s,x) 

The  change  in  responcCj,  Y(s,x),  caused  by  the  increment  Ax  is 
given  by  ^Y(ssx)  -  R(s)  AT(s,x)  (4-54) 

Substituting  Equation  (4-53)  for  AT(s9x),  Equation  (4-54)  becomes 

AY(s,x)  *  (~-\     R(s)T(s,x)S(s,x)       (4-55) 

Equation  (4-55)  indicates  that  AY(Sj,x)9  the  increment  in  the  response, 
may  be  obtained  by  calculating  the  response  due  to  an  equivalent  input 
perturbation  signal,  AR(s),  applied  to  the  original  transfer  function, 
T(s,x),  where 

R(s)  S(s,x)        (4,56) 


(-*) 


Simulation:   Equation  (4-55)  may  be  represented  in  block  diagram  form 
as  drawn  in  Figure  4-3 ,  where  S(s9x)  has  been  represented  according 
to  (4-52) «  From  Equation  (4-44) ,  polynomials  Q(s9x)  and  Q(s,o)  can 
be  written  in  exj 


Q(s,x)  ■  (c  +  xd  )  +  (Cj  +  xd  )s  +  — — — — 

+  (c  ,  +  xd  ,)sm°  +  s  (4-57) 

m°l     m-1 

Q(s,o)  ■  c  +  c  s  +  <=---  —  —  +  c  ,s    +  s  (4-58) 

o    1  m-1 
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Simulation  of  the  transfer  function  Q(s,o)  /'Q(s,x)  may  be  accomplished 
on  an  analogue  computer  as  illustrated  in  figure  4-4,  A  similar 
simulation  is  possible  for  the  transfer  function  P(s,o)/P(s,x)  so  that 
interconnection  according  to  Figure  4-4  gives  an  analogue  computer 
method  for  obtaining  the  change  in  the  response,  AY(t,x),  directly  as 
a  time  function  trace  for  any  input  signal. 
Sinusoidal  Steady  State. 

In  the  sinusoidal  steady  state9  s  s  jfe), 

T(;Ax)  -   T(U,x>  ®j9^'x)       (4-59) 
where, 

tan"1  3  -   ImT<3(,)»x?  (4-6(n 

t&ti       *  ReT(jA,x>  ^  W} 

and 

£cT(jtf,x)  -A,|T(fi>,x)|  +  je(fc)(x)       (4-61) 

From  Equation  (4-61)  and  the  definition  of  S(j£),x)  as  given  by 

Equation  (4-49)  it  follows  that; 

R«tS<».x)l  "  *t£iQ"X)  <4"62) 

Using  Equations  (4-62)  and  (4-63),  if  then  x  changes  by  an 
increment,  Ax,  the  corresponding  incremental  changes,  A^t|T(£>,x)f  , 
andA0(&),x),  are  given  by 

A&tlTfo).*)l 


-7-r—nr-    -  Re[S(jO,x)]        (4-64) 
-^SL   -  I»[S(jV)1  (4-65) 
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IV-4  -  Sensitivity  integrals 

The  expression  for  sensitivity  leads  to  functions  of  parameter 
x  (or  k)  and  s9   t  or  U>  „  From  an  engineering  point  of  view  it  would 
be  convenient  to  have  an  overall  measure  of  sensitivity  which  is 
independent  of  either  time  @r  frequency.  Such  a  figure  would  be  use- 
ful for  comparing  the  relative  sensitivities  of  different  parameters 
in  a  system  as  well  as  presenting  the  possibility  of  including  a 
sensitivity  criteria  for  purposes  of  design  or  optimization. 
(a)  -  ffhase  Integral, 

Consider  first  the  area  ©f  the  equivalent  input  perturbation 
signal,  Ar(t,x)»  relative  to  the  area  of  the  input  signal,  r(t). 
The  total  area  of  the  input  signal  is  given  by  the  integral, 


3   r(t)dt  ■  liaa   f  *  p  R(S)  ~j 
o  t*»  *■   L  s    J 


I   «■  lim 

t-><*>  o 

"   llM   „SR£s}_   .   R(0)  (4.66) 

8  ™^  O     S 

In  the  above  calculation  use  is  made  of  the  fact  that, 


J   f(t)dt  ■  AT1  I  F(s)~] 


and  the  final  value  theorem  is  applied.  Similarly  the  total  value  of 
the  area  of  the  equivalent  inpust  perturbation  signal  is  given  by 

^1   s  lim  J     &  r 
1     t-»»  o 


lim  f  A  x  \  / 
t-*oo  \^  x  /  Q** 

(A]LA      llm   %(s)s(s8X) 
>•  '       s  ->  0 

("^fv  R(©)S(o,x)  (4-68) 


LOO 


If  I  4   0,  then  the  ratio  of  Equation  (4-67)  to  Equation  (4-66) 

denoted  by  I(x)  is, 

41       /  a    \ 
I(x)  -  — j*-   -  (^f-J  S(o,x)   (4-68) 

When  I.  ■  0  it  means  R(0)  »  0  and  therefore  from  Equation  (4-67) 
it  follows  that,  hi     ■  0,  provided  3(©,x)  4  fl&.  When  R(0)  -  0  and 
S(o,x)  =  00,  then  Al,  may  b®  evaluated  by  calculating  the  limit  of 
the  product  R(s)  S(s,x)  as  indicated  by  Equation  (4-67). 

Alternately  I(x)  may  be  obtained  by  considering  the  area  of  the 
change  in  the  response  due  to  the  increment  Ax.  Prior  to  the  para- 
meter change  the  area  of  the  response  is  given  by 


I     - 

o 

lim      \      y(i 
t-»-<»    o 

t,x)dt 

a 

lim        C1 
t-»co    G* 

R(s)  T(s,x) 

8 

St 

lim 

s     R(s)  T(s,x) 

s-»o 

s 

»  R(o)  T(o,x)  (4-69) 

After  the  increment  x,  the  area  of  the  change  in  the  response  is 
with  the  aid  of  Equation  (4-54)  given  by 

t 

41   «  lim   C  Ay{t,x)dt 
°     t-»oo  Jo 

m     (A*S    lim    71     R(s)T(s,x)S< 

s  (~^~)  lim 

"  Vx")  R(0)  T<°»x>  s<°»*>  (4-70) 

If  I  4   0,  the  ratio  of  Equation  (4-70)  to  Equation  (4=69)  exists  and 
o 

is  given  by 


lim    -    R(s)  T(s,x)  S(s,x) 

8  ->0 
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(5-) 


I(x)  -  — —   -  (-a=-j  S(osx)  (4-71) 


When  I  =  0  it  follows  from  Equation  (4-80)  that  1  s  0,  provided 
S(o,x)  4   0.  When  1  =  0  and  S(@j,x)  &co  then  &1     may  be  evaluated  by 
calculating  the  limit  of  the  product  R(s)  T($px)  S(s,x)  as  shown  by 
Equation  (4-70). 

Now  consider  sinusoidal  steady  state  and  the  area  of  the 
incremental  phase  characteristic  XQ(xB&))  plotted  on  a  logarithmic 
frequency  scale s   defined  by  the  integral 

A  ©(k>9x)d(«*to) 


j  &  Qtf.%   v\    -d-^- 


A  @(Cj8x)  -SW-  (4-72) 

as 


Substituting  Equation  (4-65)  into  Equation  (4=72)  we  get 
I0(O,x)  =  J  (&■)     Im[S(.l4),x)J  4^  - 

(^/i»[s<j»,x)j 


du) 
(0 


(4-73) 


From  the  reactance  integral  theorem  (Appendix  If),  it  is  known  that 
if  H(j<0)  -  R(w)  +  jx( 

and  1  im 


then  f*  ««  d"         I 

Applying  Equation  (4-74)  to  Equation  (4-73)  we  get 

I0(x,u)  -  "K^  lim   S(s,x) 
9  2  V  x/  §^0 


(4-74) 


(4-75) 


(-4) 


Equations  (4-68)  ,  (4-71),  and  (4-75)  enable  the  following  theorem 

to  be  stated. 

Theorem  -  For  an  incremental  parameter  change,  ix,  the  product  of 

times  the  logarithmic  sensitivity  function  evaluated  at  S  ■  0 
is  equal  to 

(a)  -  The  ratio  of  the  area  of  the  equivalent  input  perturbation 

signal  to  the  area  of  the  input  signal,  provided  the 
latter  area  is  n©n°zero. 

(b)  -  The  fractional  change  in  the  system's  response  to  any 

input  signal,  provided  area  of  the  response  is  non-zero. 

2 

(c)  -  -  =   times  the  change  in  area  of  the  sinusoidal  steady 

state  phase  characteristic  plotted  on  a  logarithmic 

frequency  scale. 
When  S(o,x)  »  0^  Equation  (4-75)  indicates  that  the  increment  Ax, 
causes  zero  net  change  in  the  area  of  the  sinusoidal  steady  state 
phase  characteristic  plotted  versus  JLtQ.     Similarly,  when  S(o,x)  »  0, 
Equation  (4-70)  indicates  that  Ax  provides  zero  net  change  in  the 
area  of  any  time  response  characteristic  which  initially  has  finite 
none-zero  area. 

If  the  area  of  response  is  zero  and  S(o,x)  is  finite  there  will 
also  be  zero  net  change  in  area  produced  by^x.  S(o9x)  may  be  cal- 
culated directly  from  Equation  (4»52)„ 


c 
o 


c  +  xd 
o     o 


(4-76b) 


where,     Q(o,x)  s  c  +  xd 


103 


P(o,x)  -  aQ  +  xb  (4-76b) 

From  Equation  (4-76b)  it  is  obvious  that  S(o,x)  will  be  identically 

zero  when  ad   ■  c  b  . 
o  o      o  o 

(b)  -  The  Magnitude  Integral. 

It  is  also  of  interest  to  investigate  the  area  of  the  change  in 

the  sinosoidai  steady  state  magnitude  characteristic  as  given  by 

Equation  (4-64).  Consider  the  integral, 


; 


'oo 
AIj^x)  -  )     A^T(u,x)  d<o 


<&l 


'(A 

Re  [S(jo,x)]d(0  (4-77) 


o 
From  the  resistance  integral  theorem  (Appendix  IV)  it  is  shown  that  if 

H(j6>)  »  R(jU)  +  jx(fc)  (4-78) 

and  lim    H(jfci)  =  0 

OH00 

Then  f* 

J  R(«o)  d<o  .  |  lim    sH(s)  (4-79) 

O  s-»oo 

Substituting  Equation  (4-79)  into  Equation  (4-77)  we  get 


Ay^.x)  -  (-*f-)f   li»    «S(8,x) 


The  use  of  this  integral  is  illustrated  in  the  following  discussion, 
(c)  -  Sensitivity  Squared  Integral. 

In  many  cases  S(o,x)  is  zero  so  it  is  not  always  a  convenient 
criterion  for  comparing  the  sensitivities  of  different  parameters  in 
a  system.  This  suggests  the  possibility  of  investigating  the  areas 
of  the  incremental  error  signals  squared. 
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In  previous  sections  of  this  chapter  it  was  shown  how  the  oritur* 
ion  of  minimum  square  error  to  the  incremental  response  was  utilized. 
A  simpler  and  less  restrictive  approach  is  to  investigate  the  area 
of  the  equivalent  input  perturbation  signal  squared.  Thus,  the 
following  integral  is  deft 


2        (°°        2       f°°        2 
I  (x)  s  J    [*r(t,x>]  dt  -  J  [Ar(t,x)]  dt     (4-81) 

o  =  co 

where  it  is  assumed  that  Ar(tsx)  ■  0  for  t<  0.  Applying  Parceval's 
theoretic  to  Equation  (4=81)  s  with  the  aid  of  Equation  (4-56)  we  get 


H'  Z 

\         '  -Ice 


XS(X>  "  (M)    J   [R(«)S(s,x>]  [R(-s)S(-s,x)]ds    (4-82) 


when  R(s)  ■  I,  r(t)  »  o"(t),  then  Equation  (4-82)  becomes 

2 


J  4  CO 
S(s8x)S(-s,x)ds  (4-83) 


Integrals  of  this  type  are  tab lusted  in  the  literature  (reference  2) 

and  form  the  basis  of  analytical  design  procedures.  In  particular  this 
integral  can  form  a  basis  of  comparison  for  the  effects  of  different 
parameter  variations  upon  system  performance.  The  approach  is  ill- 
ustrated by  a  simple  example. 
Example: 

Consider  a  second  order  positional  servomechanism  subjected  to 
a  load  disturbance  torque  T  (t)  applied  to  the  torque  shaft.  The 

JL 

transfer  function  relating  th*  output  angle  to  the  .applied  torque  is 

given  by 

©  (s) 

(4-84) 


LiT< 


Where  J,  F  and  K  are  moment  of  inertia.,  viscous  damping  and  gain 


itr 


Y(s,x) 


AY(s.x) 


^x/x) 


Fig.  4-3:  Block  Diagram  for  Generating  the  Incremental  Response 
Caused  by  the  Parameter  Change,  Ax. 
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respectively.  It  is  desired  to  compare  the  sensitivity  of  the 

Impulsive  response  when  T  (t)  ■  S*(t)  for  incremental  variations 

L 

of  F  and  K.  From  the  polynomials  of  Equation  (4-84),  it  follows: 

x  ■  F 
P(s,F)  -  1 

P(s,o)  -  1 

2 
Q(s,F)  -  Js  +  Fs  +  K 

2 
Q(s,o)  -  Js  +  K 

Substituting  into  Equation  (4-52),  we  get 

S(s,F)  -  5—^5 (4_85) 

Js  +  Fs  +  K 

Similarly  it  can  be  shown  that 

S(s,K)  -       -K (4-86) 

Js2  +  Fs  +  K 

Substituting  Equation  (4-85)  and  Equation  (4-83)  we  get 


T 


and  similarly 


where 


V/  '   -j~  I  Js  +  Fs  +  K      Js  -  Fs  +  K 


4 


W    -  y/  K/J      and 

I 


2/KJ 

CO  is  the  undamped  resonant  frequency  and  3  Is  tne  damping  ratio. 
Equation  (4-87)  and  (4-88)  show  that  both  sensitivity  squared  inte- 
grals are  proportional  to  o  .  For  1t>  -r  a  fractional  variation  of 

n       v      z 

F  is  seen  to  have  a  greater  effect  upon  its  sensitivity  squared 
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integral  than  the  identical  fractional  variation  of  K  has  upon  its 
integral 

IV-5  -  Signal  Flow  Graph  and  Sensitivity 

The  general  formulation  of  the  transmittance  of  a  signal  flow 
graph  described  by  Mason,  can  be  extended  to  evaluate  the  sensitivity 
function.  Mason  has  derived  a  simple  and  very  useful  rule  for  eval- 
uating system  transfer  functions  (or  sometimes  called  transmittances) 
from  signal  flow  graph „  The  transmission  from  any  independent  node 
x,  called  sourc@pto  a  dependent  node  y9  called  sink,  is 

2p   A 

T      .   *y   *y 

xy         A 

Where  the  system  determinent  A  is 

A  =  1  •  2  L.  +  2   L°  L-;   -  2  V!   L{'  LT  +  --«•-  (4-90) 

i  l   i,j   *  J    i.J.k  *  J  k 

In  Equation  (4-90)  L  represents  a  loop  transmission,  i.e.  the 
transmission  around  any  closed-loop  and  *    L.  is  the  sum  of  all  loop 
transmissions .  L°  L°   is  the  product  of  the  loop  transmissions  of 
ny  two  loops  that  do  not  have  a  node  or  branch  in  common,  and 

^»  L'  L'   is  therefore  the  sum  of  all  such  products.  VI  VI   L".   is 

i,j     j  J 

the  product  of  the  loop  transmissions  of  any  three  loops  that  do  not 

have  any  branches  or  nodes  in  common,, 

The  definiton  of  P   A   in  Equation  (4-89)  is  as  follows: 

xy  xy 

P   is  any  direct  transmission  from  the  input  x  to  the  output  y.  Each 

A   is  associated  with  a  particular  3?  .  The  value  of  A   is  the  same 
xy  xy  xy 

af  A  except  f  r  the  removal  of  all  terms  containing  any  branches  or 
nodes  that  appear  in  P 
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Example; 

Find  A  ,  P1q,  A^  and  *   In  Figure  4=5 .  The  sum  of  the  indivi- 
dual loop  transmissions  is 

Li  *  T2  T3  T10  +  T3  T4  T5  Tll  +T3T4T5T6T12 

+  T4  T5  T6  T?  T9  +  T5  T6  Tg  (4-91) 

The  above  is  best  done  in  a  systematic  manner  by  starting  with  T. 

and  listing  all  loops  containing  T-,  then  going  on  to  T~,  and  taking 

all  loops  containing  T  that  have  not  already  been  listed,  etc. 

Next  we  find   r1.  L'  Li  which  is  the  sum  of  the  products  of  the  loop 

transmissions  taken  two  at  a  time9  omitting  all  those  that  have  any 

node  or  branch  in  common,,  To  evaluate  .**.  Li  Li,  the  list  in  *  L. 

liJ  I  j  i  i 

is  examined.  The  first  term  in  f  L.  is  TJT,T.,rt.  The  other  terms 

are  scanned  to  detect  those  that  do  not  contain  any  of  T_,  T«,  or  T... 

The  first  term  satisfying  this  condition  is  T.T-TgT-Tg,  but  Figure 

4-5  is  examined  to  check  whether  these  loops  (TJT-T.q  and  T.TJTgT-Tg) 

have  any  node  in  common.  They  dos  and  therefore  their  product  is 

not  included  in  L!  Li.  However  the  last  term  T_T,T0  has  no  branches 

i  j  5  o  o 

in  common  with  T«T3T10  and  therefore  1Jt,J£,.Jt*LS&~   appears  in 

Tl!  Li.  Next  we  note  second  term  in  L.f  i.e.,  T,,T,T_T, ,  and  examine 
L.  i  j  i        3  4  5  11 

as  above.  In  this  way  it  is  found  that, 

2  h't  V     -  <T2T3T10}  (T5T6T8)  (4°92) 

i,j 

Since  ^Li  L'  contains  only  one  terms  it  is  impossible  forZL"  L" 
to  exist.  Then 

L   -  i  -  2t  +  2  Li  Li 

4  i   i»J   i  J 

-  1  -  [T2T3T10  +  V^Vll  +  T3T4T5T6T12  +  T4T5T6T7T9 

+  T5T6Tg]  +  T2T3T10T5T6Tg  (4-93) 
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Continuing  with  evaluation  of  P.  and  **.  we  have, 

Plo  -  T^T^T^  (4-94) 

Since  each  loop  inX.L.  has  at  least  one  branch  present  in  P.  ,  then 

Aio  -  1  (4-95) 

and  PloAio  -  TlT2T3T4T5T6T7  (4-96) 

Substituting  Equation  (4-93)  and  Equation  (4-89)  we  get  T.  . 

10 

Formulation  of  Sensitivity  Function: 

Sensitivity  of  a  source  to  sink  transmit tance  T  .  with  respect 
to  a  parameter  is  defined  as 

Sx  J  -  -TO""   "    ?«/xM      <*-97) 
As  discussed  in  chapter  3,  the  sensitivity  function  may  be 
expressed  in  the  following  form: 


T 

s «  -  4 

x       F 


[» -=£=-" 


F       F' 
x       x 


(4-98) 


Where  T  (x»0)  is  the  transmit tance  from  source  i  to  sink  j  when  the 

specified  element  x  vanishes,  i.e.,  when  the  branch  x  of  the  signal 

flow  graph  is  removed  and  where  F  and  Ff  are  the  return  difference 

and  the  null  return  difference  respectively,  both  referred  to  the 

element  x.  Furthermore, 

Up   A 

X    ■  U — £i-  (4-99) 

ij         A  l    ' 

as  discussed  in  previous  sections.  Expressing  Equation  (4-99)  in 
logarithmic  form  and  differentiating  with  respect  to  &oc,  we  get 

8  *J    .    ^  pu    ij    .  iA*iA  (A.i00) 

The  terms  on  the  right  hand  side  are  themselves  sensitivity  functions 
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and  one  may  therefore  regard  the  sensitivity  function  of  Equation 

(4-97)  as  a  difference  of  simpler  sensitivity  functions.   v  .   — - 

is  the  sensitivity  of  the  graph  determinent  with  respect  to  the 

element,  x.  The  way  the  determinent  A  was  defined,  it  is  then  only 

a  function  of  the  feedback  branch,  taking  the  form  of  the  sum  of  the 

products  of  branch  transmittances.none  of  which  appear  more  than 

once  in  any  single  product.  Therefore  we  can  write  A  as  a  linear 

function  of  s, 

X*   ■  aac  +  A 

and  X  A 

$  x 

or  x  -M-  «  ax  «  A  -  4°  (4-101) 

Where  L     is  the  determinant  when  x  ■  0,  i.e.,  when  the  specified 
branch  x  of  the  signal  flow  graph  is  removed.  However  accordinging 
to  signal  flow  graph  theory, 


Dividing  Equation  (4-101)  by  A  we  get 


(4=102) 
4° 


sx  -  ixTT  ■  l  *-*'*   '  l  -  -r-  <*-103> 


Similarly,  the  first  term  of  the  right  hand  side  of  Equation  (4-100), 
which  may  be  thought  of  as  the  path  sensitivity,  is  found  to  be 

p  Z^t/ii)0  i 

x         Il*ij°ij'  x 

Then  T.  .      _      a  j       , 

s» J  -  sl  -  st  ■  -r-  -  ~h        (4-105) 

x       x 
Equation  (4-105)  was  previously  obtained  in  chapter  2.  Also  sensiti- 
vity could  be  written  as, 
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s*  ■  *  ■  z<<i  -I!)  (4'106> 

In  certain  configurations  this  may  be  an  easy  way  to  evaluate 
sensitivity  functions.  However  this  is  recommended  as  an  interest- 
ing topic  for  further  investigation  to  utilize  the  concept  of  signal 
flow  graph  and  its  related  techniques. 

IV-6  -  Multiparameter  Sensitivity 

So  far  we  discussed  the  sensitivity  of  a  transmission  function, 
T,  which  depends  upon  a  single  parameter  x  and  defined  as 

Now  consider  T  as  a  function  of  n  parameters,  i.e., 

T  -  T  (xx,  x2,  ,  xn) 

Several  definitons  have  been  proposed.  Hakimi  and  Cruz  (see  reference 
18)  define  multiparameter  sensitivity  by  , 

m*yfmiim  /.     <        \ 

o  <  a>t   <  h    H^/if,  *t)  (4-108) 

where  o  is  the  maximum  variation  of  the  ith  parameter  as  specified 
in  particular  application.  For  small  o  ,  the  fractional  change  in 
transmission  is  given  by  (neglecting  second  order  effects), 

Setting  y  °  A-x  ,  we  can  consider  the  set  of  fractional  parameter 
increments  as  a  vector, 

dy  -  [dtfcj^),  ,  d(/*xn)] 

and  Equation  (4-109)  becomes  the  scalar  product 
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At 
t 


V  [kUyv 


yn)]  *  <*y 


(4-110) 


The  sen- itivity  of  transmission  clearly  depends  upon  the  gradient 
vector  7/^T.   Kuo  and  Goldstein  (reference  18)  suggest  that  multi- 


parameter sensitivity  S  be  defined  as 
,T 


,  yn)j 


(4-111) 


This  has  the  advantage  that  it  agrees  with  the  single  parameter  defin- 

ition  in  Equation  (4-107).  The  magnitude  of  the  multiparameter  sensi- 

tivity  function  is 

|  ST  |   -   (ST*  S"T)1/2  (4-112) 

Donald  A.  Calohan  (reference  17)  defines  the  parameter  sensitivity 

vector, 

sT(ja)  m     J  hUw 


^JUsc, 


A 


(4-113) 


>i  x„     ii     J    a  %„  ii 


where 


Hence, 


i   *Xi 


oC(«*>)  a  j^|T(jw)(   ,  the  attenuation 

6  (ii)     ■  arg  T(j<^   the  phase 

x.  s  element  varying 

x  *=  a  set  of  orthogonal  unit  vectors 


5   >>* 

fc   -T X 


iXi 


Re  (8«W, 


mil 


{ 


.T(j*J) 


The  magnitude  of  the  vector  Be(S 


x 


(4-114) 
(4-115) 


),  i.e., 
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/,2  (*?-  -.) 


(4-116) 


is  therefore  an  indication  of  the  sensitivity.  Multiparameter  sensi- 
tivity could  best  be  discussed  in  terns  of  state  space,  and  is  recom- 
mended as  a  topic  for  further  work  on  this  line. 
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The  object  of  this  thesis  as  stated  in  the  "Introduction,"  was 
a  literary  research  on  the  subject  of  sensitivity  applied  to  the 
linear  feedback  control  system,  with  only  one  parameter  changing. 
Two  different  cases  were  discussed.  One  when  the  parameter  changes 
are  incremental  and  the  other  when  parameter  changes  are  moderate 
or  large.   In  each  case  a  system  sensitivity  was  defined  associated 
with  methods  of  measuring  this  sensitivity  function. 

In  the  case  of  incremental  parameter  variation,  the  concept  of 
root  sensitivity  was  presented  and  from  that  a  graphical  design  tech- 
nique discussed.  When  the  parameter  variations  were  large,  a  new 
definition  of  sensitivity  was  introduced  and  it  was  shown  that  when 
the  leakage  transmission  between  the  input  and  output  is  negligible 
the  two  definitions  are  identical.  A   design  philosophy,  and  the  use 
of  Nichols  Chart  was  also  presented.  Finally,  the  new  aspects  of 
sensitivity  such  as  specification  for  sensitivity  function  based  on 
forced  response  criterion  and  least  mean  square  error  were  explained. 
A  new  approach  based  on  the  signal  perturbation  theory  together  with 
sensitivity  integrals  were  also  introduced. 
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Appendix  I 

Bilinear  Theorem,  Return  Difference,  Null  Return  Difference 

1  -  The  Fundamental  Feedback  Flow  Graph. 

The  signal  flow  graph  for  a  linear  feedback  control  system  is 

shown  in  Figure  Al-1,  where, 

t  .  is  transmission  from  input  I  to  output  0 

t   is  transmission  from  source  S  to  output  0 

t  .  is  transmission  from  input  I  to  control  C 

t   is  transmission  from  source  S  to  control  C. 
cs 

Then  we  can  write, 

Output  0  -  t  .1  +  t  S 
r  oi     08 

and         S  -  kt  J{\   -  kt  ) 

ci       cs' 


Hence, 


0  D  t/  x  ^1*08  /A,  iX 

1  -  T<«>  -  \,i  +   1  -kt  <A1"l> 

cs 


Equation  (Al-1)  is  called  the  "fundamental  feedback  equation"  and  is 
the  most  important  equation  in  linear  feedback  theory  and  all  the 
significant  properties  of  linear  feedback  systems  are  derivable  from 
this  equation.  As  an  analytic  tool  it  shows  how  a  linear  circuit  prob- 
lem can  be  broken  up  into  two  simple  circuits. 
2  -  Return  Difference. 

The  system  transfer  function  as  derived  above  is 

T(.)  -  f  -   tol  ♦  ktclto§/(l  -  kte.) 

suppose  it  is  asked  "what  is  the  effect  of  the  feedback  around  k,  on 
the  system  function  T?"   Clearly  the  effect  is  given  by  the  quantity 
1  -  kt  .  This  quantity  1  -  kt   has  a  very  simple  physical  interpre- 

Co  Co 

tation  readily  seen  by  opening  the  closed  loop  anywhere,  for  example 
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Fig.  A 1-1:   Fundamental  Feedback 
Flow  <-rapb. 


cs 


Fig.  Al-2:   Physical  interpretation  of  loop  transmission  function, 

t 


cs 


-kt 


-kt 


OS 


Fig.  Al-3:   Derivation  of  null  return  difference. 
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at  AA'  and  measure  the  signal  returned  at  A,  the  latter  is  kt 

cs 

The  difference  between  the  injected  and  the  return  signal  is  pre- 
cisely 1  -  kt   and  is  called  the  return  difference  for  element  k. 

Fk  2  1  -  kteg  (Al-2) 

3  -  Bilinear  Theorem: 

The  transfer  function  of  any  two  port  linear  systems  can  be 
written  in  the  form  of  the  fundamental  feedback  equation  (Al-1),  such 
that  all  t   are  independent  of  the  system  parameter  k(s),  providing 
k  can  be  represented  as  a  controlled  source  S  ■  kC  and  providing  k 
appears  only  in  S  and  nowhere  else. 

The  above  theorem  is  called  bilinear  theorem  because  the  funda- 
mental feedback  equation  may  be  considered  as  a  bilinear  transforma- 
tion of  k. 

4  -  Null  Return  Difference,  Figure  Al-2. 

The  null  return  difference  is  the  return  difference  evaluated 
under  the  condition  that  the  input  is  adjusted  to  give  zero  output. 

In  terms  of  the  signal  flow  diagram  the  null  return  difference 
is  determined  as  follows.  The  diagram  is  broken  at  A  and  A'  and  a 
unit  signal  is  transmitted  from  A'.  The  signal  reaching  S  is  k.  This 
is  transmitted  back  to  C  and  also  to  the  output  0.  If  the  input  is 
adjusted  to  make  the  output  zero  the  signal  -kt   must  be  arriving 

OS 

at  0  along  the  branch  from  I.  Hence  I  is  -kt  /t  ..  This  input 

os  oi 

is  also  transmitted  along  the  branch  from  I  to  C  and  the  return  signal 
is  the  sum  of  the  two  signals  arriving  at  C.  Then  the  null  return 
difference  with  reference  to  k  is  given  by  the  relation, 

F.  •  -  1  -  (kt   -  kt  t  ./t  .) 
k         N  cs     os  ci  oi7 
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The  system  transfer  function  can  be  written  in  terms  of  return  diff- 
erence and  null  return  difference, 


kt  ,t 
T   «   t   + £i_08_ 

1      oi  T  1  -  kt 


n    £ 


cs 


oi 


kt   +  kt  .t  /t  . 
cs     ci  os  oi 

1  "  ktcs       (il-3) 


Where  P'  is  the  null  return  difference  and  F  is  the  return  diff- 
erence with  reference  to  k„ 
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Appendix  II 

The  Curvature  of  the  Root  Locus  at  Ordinary  Points 

The  curvature  of  a  curve,  which  is  defined  by  means  of  an 

implicit  relationship  F(x,y)  »  0  is  given  by 

,                      FF2-2FFF+FF2 
1  xx  y xy  x  y         yy  x 


where 


F 


(F  2  +  F  hV1 
v  x     y  ' 


^F_ 
ax 

by 


XX 


xy 


ox 


F 


yy 


*y 


0 

ex 

*y2 


V>yJ*  «■ 


F 

x5y" 


(A2-1) 

(A2-2) 
(A2-3) 
(A2-4) 

(A2-5) 
(A2-6) 


If  we  introduce  a  new  cartisian  system  with  its  origin  at  s.  and 
oriented  so  that  the  tangent  to  the  locus  coincides  with  the  u-axis. 
Then  in  this  system  for  the  implicit  function  F(u,v)  we  have 
(Figure  A2-1) 


and 


*F 
bu 

v-0 
u»0 

-     Fu(0,0)     -     0 

1 

-     - 

F     F2-2F     FF    +F     F2 

UU   V                 UV    U   V           w   u 

J> 

rt2JL,  2x3/2 
(F       +  F     ) 
u             V 

(A2-7) 


and  since  F  (0,0)  ■  0  and  curvature  is  evaluated  at  point  (0,0),  then 


1 


F 


uu 


u«0 
v«0 


The  center  of  curvature  is  on  the  line  normal  to  the  curve.  Because 
of  the  perpendicular  choice  of  coordinate  the  positive  sign  of f  means 
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S-Plane 


Fig.  A2-1:   u»  v  coordinate  system. 


Fig.  A2-2:   (S  -  S  )  in  v,  and  u  coordinate. 
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that  the  center  of  curvature  lies  on  the  positive  v-axis  and  vice 
versa.  Since  the  gain,  K(sX  is  a  rational  function  of  s  and  analytic 
except  at  its  poles  s  »  z  ,  then 

n  m 

M*)  »/n,[K(s)]  -   2M«-P.>-  S  JM«-Z4)+t\J 

J-1      j    «        (A2-9) 
This  is  a  two  dimensional  complex  potential  function  such  as  would 
be  used  to  describe  the  field  around  a  set  of  n  unit  positive  line 
charges  at  the  points  P  and  m  unit  negative  line  charges  at  the 
points  Z  for  such  configurations,  the  lines 

Re  K  (s)  =^|k(s)|  ■  constant  (A2-10) 

would  be  equipotentials,  and  the  lines 

Im  K  (s)  -  argK(s)  -  constant  (A2-11) 

would  be  stream-lines  or  lines  of  force.  Since  the  above  relations 

and  the  relation 

m  n 

arg  K(s)  -   2  «*8  («  -  Z.)  -  £  arg  (s  -  P  ) 

i-1  1         j-1  X  (A2-12) 

are  equivalent  with  a  particular  value  for  the  constant,  it  follows 

that  the  root  loci  coincide  with  certain  lines  of  force  in  potential 

analogy. 

Having  the  above  idea  in  mind  we  can  define  function  F  as  F  • 

ImK  and  sometimes  is  defined  as  F  ■  ImoitK,  which  is  more  suitable 

for  our  purpose.  As  K  takes  on  positive  real  values,  F  is  equal  to 

zero  and  therefore  serves  as  the  real  function  of  Equation  (A2-1). 

Therefore 


(A2-13) 


,                               F 

1                      uu 

"35"                   F 

<s                           V 

Im 

(ImAc) 

'uu 

r        1 

Im 

'<MDV  1 
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In  general  K  is  expressed  as  a  function  of  the  variable  Z  ■  u  +  jv 
and  since  InK  is  an  analytic  function  then, 


»u    " 

tr 

-fe-- 

.  \          d 

3  \jv             J  dZ 

Hence 

1 

Im(^K) 

Im(^K) 

J>   " 

Imj(^K)z 

Re(ifJC)z 

To  evaluate  Equation 

(A2-15)  we  assume  Z 

■  (s  •  s  )e  " 
o 

function  f  we  have 

df 
dZ 

df     ds      df 
ds     dZ      ds 

ejy 

d2f 

2    /   \2 
d*f   /  ds  V 

df    d2s 

dZ2   *   ds2   VdV   ' 

d«   dz2 

d2f    2J* 

■  -   0   e 

ds 

and  similarly 

dnf 

dzn 

dsn  V. iz ) 

J   n 

ds 

(A2-14) 


(A2-15) 


tt 


(A2-16) 


ejn*     (A2-17) 
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Appendix  III 
Derivation  of  an  Expression  for  the  Sensitivity  of  a  Multiple  Order 

System  Root 
If  we  write  the  expansion  of  the  total  differential  of  L(s)  to 
include  higher-order  terms  (see  reference  15,  page  81),  we  have 


(«.), 


where 


.  o  -  ^ds  +  4£-dK+i^dz<  +  r^-dpJ 

1  8^  *■       oK       i  *Zi   l   i  *Pi   *•  I 

2,   I  same  as  above  bracket! 

♦tK  [ ]      *  

Ju-  (di)-l«  +  *i|^l  J£—   (d.)-2(dK)2 


+  .  --* 


^s*~hvL  2'      ls*~hK 


+  -ii 


(dK)a 
OK 

Next  retain  only  the  lowest  order  terms  for  each  parameter  and 

note  that  the  first  (N-l)  derivatives  of  L  with  respect  to  S  are  zero 


at  S  ■  S   .     Then  , 


•'!    •      [<"l>m  "]  [l 


*■    »  +  Z  -ii-  «.  +  J  YL 


ir  * + f  nr  dzi + 1  ti^-  dpi 


This  suggests  the  following  notation: 

r  •  i m 

dsi  -  [* -T- +  *  sz;  *. + j  *P; dpij 

and  comparing  the  above  two  equation  we  have 


1/N 


J  •-• 
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K 


(-DN  n: 


8    ■    8 
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Appendix  IV 
Resistance  and  Reactance  Integrals 
To  prove  the  resistance  and  reactance  integrals,  we  start  with 
Cauchy's  residue  theorem.  Consider  any  region  defined  by  a  closed 
boundary  in  which  and  on  whose  boundary  H(s)  is  single  valued  and 
analytic,  except  for  poles  of  any  finite  numbers  of  multiplicity. 
Then  the  line  integral  of  H(s)  around  this  boundary  is  related  to 
the  residues  of  those  poles  of  H(s)  which  are  located  inside  the 
region  as  follows, 

&    H(s)ds  -  2Tij  ^ residues. 
If  H(s)  has  no  singularities  in  the  right  half  plane  and  as  on 
Figure  AA-1,  then 


<£    H(s)  ds  -  11m     J    H(s)d(ju»  +  lim  J   H(R,9)d(Re"J0) 
J  R-»»   -jR  R-*»  71/2 

-  lx    +  I2  -  0  (A4-1) 

Next  we  consider  some  properties  of  the  loop  transmission  function 
H(s).  The  loop  transmission  function  can  be  regarded  as  transforms 


of  impulse  responses  of  real  systems.  The  transform  of  a  time  function 

»8t 


is 


H(s)  -  J  h(t)e"< 


*dt  (A4-2) 


-oo 


So 


J     h(t)   [cos&Jt  -  j   sinOt]   dt 


where 


H(j6»  - 

-oo 

-  R(w)  +  jx(w)  (A4-3) 

R(fc))  *  J  h(t)  cos^Ot  dt  is  an  even  function  of  C*> 

-oo 

and         x(W)  "°°  J    h(t)  sin  CO  t  dt  is  an  odd  function  of  d>  . 

-oo 
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S-Plane 


Fig.  A4-1:  Application  of  Cauchy's  residue  theorem. 


S-Plaoe 


Fig.  AA-2:   Contour  for  F(s)/s. 
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Thus,   H(jto)  -  R(tJ)  -I-  Jx(6J)  -  H(-j«J>  -  A(-U))   -  Jx(-<0)  (A4-4) 

It  is  often  necessary  to  focus  attention  on  R(s)   for  s  near  zero  and 
for  s  near  infinity.     Expressing  H(s)  about  the  origin  gives 

H(s)  -  A    +  B'    s  +  A"  82  +  B"  s3  + —  (A4-5) 

x/oooo  N/ 

where  all  the  A'  's  and  B'  ' s  are  real.  The  notation  for  the  expansion 

of  H(s)  at  Infinity  is 

H(s)  -  A^+  (B^/s)  +  <AV*2)  +  <A4-6) 

where  all  Aj^  and  Boo  are  real. 

Consider  H(s)  so  that  it  satisfies  Equation  (A4-4),  is  finite 
(or  zero)  at  infinity  and  is  analytic  in  the  right  half  plane.  Apply- 
ing Cauchy's  theorem  to  [H(s)  -  H  (»)]  over  the  right  half  plane 
whose  boundary  is  the  ju-axis  and  the  right  half-infinite  semi-circle, 

£   [H(s)  -  H((0)]ds  -  0 

J  (R 

-  li«  ;   [H(j*J)  -  AJ  jd^> 

R-»o>  -R 

+  lim  J    [H(RejQ)  -  AJ  d(ReJ0) 
R-w»  -Tl/2 

(A4-7) 

But  let      H(jw)  =  R(u>)  +  jx(<0)  ,  and  note  that, 

R(»)  »  R(-63) 

x(tf)  -  -x<-6>) 
Consequently, 

lim   j   [H(jU)  -  AJ  jdco  -  j2  j  [R(«J)  -  AJ  d(^ 
R-*»  -R  0 

As  8  approaches  infinity  the  only  term  of  [H(s)  -  ^J  that  contributes 

to  the  second  integral  in  (A4-7)  is  (B^/s)  and  then 

(*/2  -TJ/2 

lim   /     [H(s)  -  AJds  -  lim  j     B^  /s  ds 
r-*co  -n/2  R-*»T/2 
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CW2  19 

-   lim   )  B'^/s  ReJW  d9 

R->*  7T/2 

-J/2 
R*oo  T?/2 


and  therefore, 

J2  J  | 


[R(60)  -  A  ]  d<0  -  jlTBjo  »  0 
0 

or  .» 

J  [R(«)  *  A  J  dU>    »   VaiB^  (A4-8) 

o 

Note,  that  if  we  had  used  H(s)  rather  than  £H(s)  -  AJ ,  as  the  integrand, 

then  the  integral  will  be  infinite  unless  A^-  0.   If 

11m  H(s)  »  A_  =  0 

S  ■»«© 

then        B1   -lim  SH(s) 

S-»  oo 


and        J 


o© 

R(W)  d60  -  l/2TTlim  SH(s)  (AA-9) 

0  s->  oo 


This  is  the  simple  form  of  "Resistance  Integral"  used  in  Chapter  IV. 
Reactance  Integral: 

We  apply  Cauchy's  theorem  to  H(s)/S  over  the  contour  Figure  A4-2 

-TT/2    A 

ir/2 


£  afai  d8   .   lln  (7  jll^l  +  /"*  J^L  d(ReJe) 

J     8         t*o»      U       -ff/2     ReJ0 


R-*co\ 


+  J'   JliM-  do)+  J7'   -4V  «"jf>(  (A4-10) 
<R     <*>         <TT72    reJ^         J 

When  the  first  and  the  third  term  of  the  above  equation  are  combined, 

the  odd  parts  cancel  and  the  even  parts  add.  Writing  H(j«>)  -  R("J> 

+  J*(W),  Equation  (A4-10)  reduces  to 
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;aO 
x<w)  -&g    +  n  <ao  -  AJ    -    0 

—  CO 


■  CO 


J  x(W)  -^       .     „fl?2   (Ao  -  AJ  (A4-11) 


-co 


If  A      -     0,   then  Equation  (A4-11)  reduces  to 


J     x(W)  -~^      -       -17/2  lim       H(s) 


-CO  W  8*0 

which  is  the  simple  form  of  reactance  integral  used  in  Chapter  IV, 
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This  thesis  is  the  result  of  literary  research  on  the  subject  of 
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design  application. 

(2)  Definition  of  sensitivity  function  applied  to  the  large 
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